ED 1U3 511 

I 

A'DTHOR 
TITLE * 



INSXITOTIOH 

SPONS AGENCY 
PDB DATE ' 
NOTE 



EDRS PRICE • 
DESCRIPTORS 



IDENTIFIERS 



DOCOMENT 'BESOflE _ ' < ' ' 

' . * \. , _ \ °22 988 " 

Charp/ Sylvia ; And ethers' / 
Algorithms^ Computation, and Mathematics (Fortran, 
Supplement). Teacher^s Commentary. Revised . 
Editioji., ° , . 

Sta-nfcrd Oniv., Calif. School Mathematics StUfJy * ' 

Group. ' , ' ' r 

National Scien'Ce Foundationv Washin*gton, D.C. 

104p.; for related -documents^ see .SE 022 98:5-987; • N<y^* 
available in hard Copy due to marginal legibility of T 
original document i • • * . . . ' - ^ 

« ^ ' ^ 

. llF-$0,83 Plu^ Postage. HC Not Available .from JIDES. 
^ *AlgorithiDs; ^Computers; *Mathematics Education; 
Programing Languages; Secondary Education; *Secondarx 
Sqhool Mathematics; -^Teaching Guidesj 
4 ^*FORTEAN; *School Ma;thematics Study Group 



ABSTRACT • 

^ • This is the teacher's guide and commentary for the 

SMSG^ textbook Algorithms^ Computation, and Math^ematics ([Fortran 
Supplement). The 'teacher^'s commentary provide^ background ii^formation 
'for the teacher^ suggestions for activities found in the'^Fortran 
Supplement, and'^answer^, f or .exergis.^:^ and activities.. The paursa is 
designed for high school students in grades 11 and 12. Access to a/ 
computer is highly recommended. (*RH)* 



* ^ Documents -acquired by ERIC include .Hany informal, unpublished * 
materials not available^ from othei^^spiirc^s. ERIC makes every effort * 

* to obtain ^the best* copy aVailable^^-^Nevextheless, items ofmarglnal * 

* reproducibility are' often 'encountered and this affects tt^ quality ^ 
of the microfiche and hardcopy te productions; ERIC makes available * 

^ ^ia cUre 'ERi'C Decumbent Reprodiictidn Service (ECfiS) . EQRS is not ' , W 

* responsible for t6.e qua'^li-^y of the. original document. Reproductions. * 

* supplied -by EDRS are the* best that can* be m^^de from the original. * 

jjc * j|6 :#t :flc :««<e :^ ******* :jc * * i * 



c. 



ED 143 511' 

A'DTHOR 
TITLE 



INSTITOTIOR 

SPONS AGENCY 
PUB DATE " 
VNOTE 



EDRS PRICE • 
DESCRIPTORS 



IDENTIFIERS 



DOCOMENT 'BfiSOUE 



SE 022 988 ' 

Charpr' Sylvia ; And Others' ^ 
Algorithms^ Computation, and Mathematics (Fortran, 
Supplement) • Teacher»s Commentary, Revised . 
Editio;i», 

Sta-nfcrd Oniv., Calif. School Mathematics StUfJy ' ' 

Group. , , 

National ScieirCe Foundationv Washington, D.C. 
66 " ^ ' S ^ " , 

104p.; for related -documents*, see J5E 022 98:5-987; • N<y^' 
available in hard Copy due to marginal legibility of \ 
original docu.meB,tv ' • ' . , . ' . ^ 

^ * * 

llF--$0.83 Plu^ Postage, HC Not Available .from JEDES, 
♦Algorithms; *Computers; *Mathematics Education; 

''Programing Languages; Secondary Education; *Secondar]t ' 
Sqhool Mathematics; -^Teaching Guidesj 

^♦FORTRAN; *School Ma;thematics Study Group 



ABSTRACT - — 

^ • This is the teacher's guide and commentary for the 

SMSG textbook Algorithms, Computation, and Math^ematics ([Fortran 
Supplement}. The 'teacher-' s commentary provide^ background ii^formation 
■^for the teacher, suggestions for activities found in the'^Fortran 
Supplement, and answers, f<^r .exergis^s^ and activities.. The pourse- is 
designed for high school students in -grades 11 and 12. Access to a . ' 
computer is highly recommended, (*RH)' 



ejIc 



:4e :te:4e :4e * :4e * :4e :4e :4e :4e :4e :4e :4e :4e :4e :4e :4e 4c :4e :4e :4e ♦ 

* ^ Documents -acquired by ERIC include niaHy informal. unpubMshe 
f materials not available from othep^SL04irc^s» ERIC makes every ef 

* to obtain ^the best* copy ayailable,^ .Nevertheless, items of margl 

* reproducibility a re* often | encountered and this affects tl^ gual 
of the microfiche and hardcopy teEroductions; ERIC makes availab 

* Via 'ERfC Decumbent Reprpdiictidn Service (ECfiS) . EQRS is not ' 

* responsible for t6e qua^Li"^jf of the, original document* Reproduct 

* supplied -by EDRS .are the- l^^st that can' be m^^de from the <2>rigina 

. • ' • 'A- 



d * 
fort 

nal * 

ity t 

le * 

ions * 

1. ' * 

^ sfc ^ jfc ^ 




7 - 




ALGPRITHMS, ^ 
COMPUTATION^^'' \ 
'AND "V / 

Mathematics . ' • 

(Fortran Sup{)!empnt) 



Teachef s Commentary \ / ' 
Revised EditiofL 



■/' ^ 

The following is '^list^f all those who 'participated in the preparatfon of this 
volume. * ■ I ' . * ^ ' ; 

Sylvia Charp,^nbbbins Teqhnical High School/Philadelphia>enJ^lvnnia 
. Alejjandra Fohythe. Gunn Higfi School, Palo A(to, California ' I > , * 
B&rhard A. G^41er. tJniven>ityof Michigan, Ann Arbor, Michigan ' ' \ 
John G^crriot, Stanford Ui;iiversity, California 

Thonias E. Hull, Uniyersity of Toronto, Toronto, Ontario, Canada 
Tho;nas A. Keenan, University of Rochester, Rochester, New York 
Robert E. Monroe, W^yneState Uni^rsity, Detroit, \Udiigan 
-5ilvio-0. Nayarro, Unfjgersity pf-Kentucky, Uxington, K%tucky " 
^ EllfeuJOrganick;University9fH<Juston. Houston, lexas ^ 
' Jes$e Pectui^iham.^^^ Unified School rXi'Mrrct, Oakland, CaliCornia 

George A. RoUimon, Afgonne National Laboratory, Argonne, Illinois' 
Phillip M. ShermaLi/Bell Telephone Laboratories) Murray Hill, New fersey 
Robert E Smith. ControKData Corporation, St. Paul. Minnesota ' ' . 
Warren Stenberg, University of Minnesota! Minneapolis, Minnesota ' ' 
Harley TilUtt, Ut S. Naval Ordnance "lest Station, China Lake, CalifSrnia 
^ . LyneyewWaldrop, Newton Sou^h High Sriiool, 'Newton. Massachusetts 
The follQwing were the principal ronsultaiVts;; ^ ' 

George E. Forsythe, Stanford University, California * ^ 

-"Bernard A. Galler, University of Michig:m, Ann Arbor, Michigan 
Wallace Givens, Argonne Nationa^l Laboratory, Argonne, Illinois - > 



■ 3 



/ < 



> 1965 an<tl966 by.The BoAtd of Trustees of the Lcland Stanford Junior Univeraity 
' ' \ All rights reserved . 

Printed i\ the United States of'Araerica 7 




Permission to make verbatim use of matermin this hook tfiust be secured 
from the Director of SMSG. Such permission will be granted except in 
unusual circumstances. Publications incorporatingSMSG materials must 
include both ath ackriowledgmexitoj tlu SMh}^cppyright^(Y^le Unii^ex- 
sity or Stanford University, as the case may befand^a disclaimer of^SMSG ^ 
endorsement. E^cclusive license will not he grmted save in exceptional 
circumstances, and then only by specific action i^f the AdvisoryBoard of 
SMSG, . . 



( 




Financial support fot the School Mathemaiics Sm4y\y(0up has been 
provided by ^he'National Science Foundation. | I ^ 



A' 



TABLET .OP CONTENT^' 



"J 



Chapter* 



IMPUT-OUTPUT MD* ASSIGNMENT ^'TATfif^NTS 
' SuHimary of Chapter P2 

F2-1. Introduction.. ^ 

F2-2. - The Character Se^''-'. , . />. . . 
_ ' ''Answers to'i^c^clsefe*.* 

.F2-3. Input-Outpirh S.-eat^ments> • ^ 

Answers to Extercls'es . 

Answers to 'Exercises F2-i+ *. . 

Ar>6wers,to Exercises F2-6" *. ' * 

Ans^^rs to Exercises 52-7*. 

. BRANCHING* A^flrSimSCRira \f\ ' • \ 
••Introduction. . • 

Outline Chapter F3 ... .\ 1 

Answers to. Exercises F3-1 . . 

Answers to Exercises F3-2 ' . . . ^. •. . \ 

'Answers , to Exercises F3-3 4 \ 

P3-^- ^P^cedence Levels for Re^atioris,. . 

'Answers t€> Exercises F3-5 ,* ^ . . v. . *. 
' ^^swers to^&cercises .F3-6 . 1 «... . 




LOOPING 
• Summary of Chapter pl+' . . . 
/Answers' to' Exercises pi^-l . 

Answerjs to Exercises Y\~2 ^. 

An'^wers Exercises Fl+^- j- . 

An^w^rs to Exercises Fj^-1+ . 



SUBPROGRAMS^ • • \ 
Use: of the Computed -00 TO 'Statement 
•Answers td^Exercises P5-I 

Answers to ExerciseS F5t3. 

Answers to Exercises , , 

^ Supf>lementary* Exercises for Section 5 -1+ 
Answers to Exercises ^P5-5' . . 
Answers to Exercises 'P5-6 - . 



19 

20 
'22 
'27 

3\ 
37 
38 
Ul 



U3 
58 



65 
66" 
'67 
73 
'79 
aif 
86 



SOME MATHEMATICAL APPLICATIONS " * ' ' 

Answers Exercises F7-1 89^ 

Answers to Exercises F7-2 ! 92 

•Answers to. Exercises i7-3 ^ . .* !• 9!^ 

Answers to Exercises F7-^ ^. 48 

Answers^ to Exercises F7-5 '.^ . * .. f . ^ ^ ! 9j 



> * Chapter ' TF2 
INPUT-OUTPUT AND ASSIGNMENT STATEMENTS' 



Summary of ' (Chapter |^ ' ^' * ^ ' 

. • ■ 'ip' . ' ^ . • ' 0— ' \ . 

Enough of the FORTRAN language, is ^ntrx)duced in this chapter to^ enable ^ 

a sifcuderit to write v^iy simple programs. One series .of exercises 4.s arranged 

in such a way as to build up complete programs-, any or all of which Xikn be . ' 

computer tested." These are Exercises 1 - 6 at the end of Sections F2-3 

Se1> A^. F2-3^ Set C, and F2-7. ' ' . 

V In addition Exercise 8/Sectl<?n PS-T, /-s also recommended^ for- computer 
.testing. It will be interesi^Aig to ^the bettei? students,, You may need to 
give the students spme specU^ h^lp 'with this one, ' ] - ' 

• . ^ ' - . ^ % . ' ' 

The outline . - * 

F2-1 »Some background on What is FORTRM ^and Vhat FORTRAN programs ^ 
look like. ' • ' - 

\ • 'F2r2 The elements bf the language, its* characters, numerals for 

^ constants ^d ■ ^^t.ement. labels, yaAables , ^ aames -for functions/- , 
and operators . An ^important^ new ,'idea to be qn the lookout for 
,^ . isl FORTRAN* s distinction, between integers arid real types ^pr 
nuiAbers and variables,- ' ■ 



\ 

F2-3 



.The! READ, PRINT" and FORMAT ^statements are int reduced T^ormat 



/ 



5-S 
of 
ind 



ntrpdticed piec.emeal at, first; only , the simplest concepts 
*05^t are discussed (l; and - F fie.lds and ret)etiti'on 
.cato3rs, Arf ields ^are' intfro^uced at the 'end of the chapter); 



Eaci^ succeeding chapter ^adds more material. 



The 
has 



assignment-restatement^ is. explained largely in terftiB tDjT^what' 
been learned Trom the extensive- ftateriar in the flow chart 



te:<fc. New ideas peculiar to^FORTI^ jmich th^;student*must be i 
onri^e lookout for are: 'the notions of prohibition 6f Mixed-' . 
. mxVe expression^, and integer division\in FORTRAN^ Vhich 



explained in terms of the^ greatest integer function. 



TF2 ' ' ' ^ 

F2-5 'Order of* ccanputation in a FORTRAN expression is, explained in 
. tenns of tlje material given in, the flowchart text.t5 

\ F2-6 . Converting integers to feals and vice versa is shown to be 
acccxnplished by. the assignment, statemeittffl 

A simple (but ^onrplete) 'FORTRAN program is displayed and a set 
of exercises given ^here the student is asked -to write his 
first complete programs/' , . * 

^ F2-8 , Presides some clerical details ifi preparing punched card for 

program statements: the continaati6n c ode^ and" the' irreleyaii^ce 
of blank spaces. *" , . ' 

. <^ * 

F2-9 ^Carriage •control for yertical spacing^ of printed^ r^siilts ,is 

^ introduced. Discussion continues in Chapter 3* This sectix)ri * 
pertains only to tl)ose Glasses who are using line printers for 
. . the output device*., With time<sharing and remote typewriter 
.consoles becoming available to an increasing extent the topic ' 
of carriage control wiH become unnecessary in a beginning course, 

F2-10 ^^lains how to read into memory alphabetic data'iin contrast to 
numeric data. The A- field format code is' xised. If. time is 
' • short, this section can be omi'tted temporarily. It is applied 

^ at "the ertd of Cha?p'£er'-5 in the exanrpiles'^^on "s^trlng manipxilation" 

* / and again in Chapter '.8. ■ ' • r 
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Literature ' £n FOBTRM II , 

- * * ' / *\ " ^ 

* , " A»non-exhaustive list is provided heie. 

A. Reference manuals , (The,- manuals Hst^d_,are revised frequently and 



the doctiment. numbers ths(t,are given may not reflect tRe latest 

revision, ) , . ./ * . * n ^ 

f ' ; * / . Computer 



X. 


"709/7090 FORTOM" 

IdM isuiietin i/cio - - 


IBM 70^ 7090 


2. 


• 

"7070 Series Programming Sjcsteme" 
im Bulletin C28 6170 ' . 


-^7070, 707^^ 


3. _ 


IM 1620 FORTRAN II specifications . ' 
im Bulletin J26 - 5602 .1 . 


^620 


h. 


"PORTRAN* f or the IM lUOl" 
im Bulletin J2h - 1U55 - 0. 


lUOl 
** 


5. 


"FORTRAN *fpr the l^flO" 
im Bulle-^in J2h - ll|68 - 0 


lUlO ^ * 


6. ■ 


"FORTRAN System for the TJontrbl Data 
160U Cpmputer" Publicatioa 087 A 


,€DC 160U 


7. <r 


"160-A PdRTRAN/Referea'ce Manual" " 
Publicqttion 60051302 


■ 160 A 


8. 


^Honeywell 80O* Algebraic* Compiler 
ManuaX" Bulletin DSl - hhk ' ^ 


Hvf*800 . X 



Primers ^ guides and other texts, 7 * 

* • • * ' " . ^ ' ^'^^ . . ' 

The l^est source of such material Is the SMSG annotated bibliography 

entitled "Study Guide in iigitai Computings and Related Mathematics^," 

which is reprinted at the end of* the Teachers Commentary for the , 

/ 

Main Text. See especially the references mentioned in Section JII, 
. Algorithpiic L^guages. . i 
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F2-1 Introduction ' ^ ' . ^ ' ' ' 
1_ , 

The FORTRAN statement* * and the declaration ' ' 

The steps of a FORTRAN program that correspond to the event or action 
boxes an a flow chart are called statements . To be complelCe, most TORTRAN 
-programs also include certain statements. which are descrip^ive«> these iiave 
'no direct counterpart in the flow chart lan^yiage. j In Section F2-3, the . ^ 
FORMAT "Statement" , is re\explained in the' new context. It is-- descriptive 
or dec;j,,arativer It is riot an event ^ or action step. At this point, we chose 
to harve the student call such nonexecutable statements as declarations rather, 
than statements. This distinction is then maintained in the remainder of 
the manual. ^ »' * ' ' 

■ ''^^ ' _ ' . • - 

" Target programs and source programs * ~ . ^ 

Terminology changes rapidly in a field which i^ moving as rapidly as 
the computer field. Most of the older literature used J:he term object 
program • Wa^r.^ using the cur rent ly* preferred term in choosing to use th^ 
word t^arget , J ' ^ . * ^ . - , - 

In additior> to the normal manner of processing a source program as 
described iij the student text, there is an alternative approach Vhich_is' 
worth knowing" a'boui^^^In this apprd>ach, "Vhe processor or compiler' program, ' 
produces ^ target 'program whiph is executed in the " interpretive~"mode I^^ 

^ This type jdf target program consists of instfuctions that are not * 
strictly machine .codes . They are machine - like instruction c^Qde^, often 
^called an V interpretive" -code. • * ' 

In order to execxite such a target a specially developed "interpreter" 

program must b^ stored in memory along with the ""target before execution can' 

proceed. 'Such compilers have been veiy successful, especially jDn machines 

• " • ^ * , f ^ 

with limited memory such as the IM l620. The interpreter program has. the 

taslc of interpreting and then carrying out the intent of each pseudo 

. instruction .of the target code. The interpret^er program in. a sense ^simulates^ 

a computer within a computer. The success of these compiled is^-explained 
^ . ' * * , -* • ■ • 

' by the fact that a FORTRAN* source program translated into interpretive code 
often occupies far less memory (fewer pseudo instructions) than a FORTRAN ' 
source program which is translated into machin^ code (mdre actual instruc- 
ti'onsy. In the former (interpretive! case, the total combined memory requiue- 
ment for^the interpretive code produced by the compile^r anji tTTe interpreter 
program is normally less, than that for the straight machine code. ' It i^ for 



• / 



this reason that this approach is popular for machines of limited memory. On 
the other hand, the latter (machine /cQde)' ca^ normally results in faster- 
running programs; • 

^ Some computers haye been desygned ariJ^built so that ^the task of compiling 
is made easy.. The Burroughs -3-55/0 computer Is an example, ,An AlLoL -compiler 
for this computer develops a target program expressed in machine code 'vhich' 
is 'as compact as could be obtaii^ed with most interpretive approaches— so the 
advantages* of both storage economy and rvmning speed are thereby achieved ' 
with essentially^ none. of the disadvantages. 

? - * • 



• / 
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F2-2 The character set ' ► » ' 

' ^ r • 

The characters shown in TalDl^ F2-X are atLso shown in the card picture,. 

Figure ;L-15 in Chapter 1 of the main'ttext. 

The special arrangemerrt of letters and jiigit^ in Table FS-l is to 
emphasize that letters and digits in the safae column of the table have one 
holer pvuich of their code lin cotranoh. For example, if you inspect the card 
picture\;^u will see that the letters E, N, and V* and the digit 5 each 
«hav6 a punch in- Row 5, . • ^ j» *^ v 

Ohara>2tef Row punches used 

E „ ' ' ; . 12, 5 

' N ' '\ ' 11, 5 ' _ ' . 

' . ' V . - ) . 'o, 5 . 



Special characters - * - . 

Many of the special characters in Table F2-1 werei placed on key punches 
when FOI^TRAN became available seve_ral years ago. Prior to -that time the 
special characters preferred by the business community wer^ used. The 
^bu^intss cothraunity was and still. is the? largest user of key puifthes.' So, ' 
i^ you need key^ punches for your aaboratory' classes and obtain the lise* of 
a "business" key punch you can expect the keys to display a differefttt set. 
of special characte^rs-i There is an equivalence between the FORTRAN' set*, 
and any other in the sense that up to how .all key punch machines of theVIBM 
Modjl 026 variety punchy only one set of hole com3binations--l.e. , the ones show 
on the card picti^e. Only the characters printed on the keys, arid^ the corres- 
ponding characters that print at the very top of. the card in.^6a^h column, ma^ 
differ. The^most typical equivalence is 

FORTRAN • Business 



+ 



) " ^ (lozenge) , 

; . , . \ - . <g (ampersand-y 

. , . ' @ (at*,*^eactty 

(apostrophe)^ (dash) • 



,Ansyers to exercises 'Set A 

l.p -.02E2 ^ : ^ . s' 

-.002E3 ' 
, -•0002Ei; ' * ' 



Assunjing 2.^h om to thS inch, thpre are ^•5V x 12 x 528o, ' or 
I..609 X 10^ cftTiiile. • , 

1.587 X 1.609 X 105 ^ •= 2.55^ X 10^5 

In "E"' notation this is 2.5'33E15 . . ■ ' 



3.. i^^6, 2 X lo5, aCs, -522.lt, 5 X 10^ 
^ . h. 3'.9lE-5, ■9.09EO/ r6.67E3, 176.!tE-9 ' 
5, None do. \ ^ - 



* 1,79 X 10 

* 6179 X 107 
_^ . 16.79 



-3 



17.9E-1+ ' or l,'Y9E-3 
.6179E-2 or 6179. E-2 
16. ^9 or 16.79EO ^ 



Statement labels ^ ' • \ * , . ^ ✓ 

^ Although 5 columns are allowed for punching a statement label, oddly 
enough, not all FORTRAN processors will accept labels as large as 99999. 
The largest^value for a label which is Acceptable ,to all proce^so»s is 
^768. If the ^ label 'has ^Less than 5 digits, it may be punched anywhere 
*in the label field, (Columns 1 -^); i,e,, tlanks are ignored,* 



Blank Characters in an Expression ' 

In F^ORTRAN II compilers spaces are ignSSjgd when the'i^ appear inside an ' 
arithflfetic expression. Hencd spaces appearing within any component ' of ^an 
e^^resslon are also ignored. Thus ABLE and ABLE are .considered . 
equivalent by the FORTRAN, canpiler/ It. happ^ens.that there are some other, 
computer languages and their processors, including some versions of FORTRAN IV 
and PL/I, where blanks are not igribredr. In* these languages .blanks ate used 
as. separators between the conrponepts of an. express ion.- 



ERIC 



sen tl 
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It^is ^therefore a* good idea to avoid using spaces "^inside variables or 
constants. In this way our rules for constructing these components are in 
harmony/ with -the other Iphguages. ^ ' , ' * 

Type of number represented by a variable rtame 

Jt might be^ interesting for students -to learn that, in more advanced' 
versions of FOMRAN, especially FORTRAN IV, a useful technique known, as 
^ -explicit typing is available to break^the leading-letter rule for deciding 
the variable type, "^e" tecliilqtie a'i^o5^r*ti5 .to declare the .variable type by 
a s o - c alled ^' "^pe : (^.c laj^aiFion . For example, ' ^ ^ . 

and REAX « ' I, j'. JJ 

. are' type declarations whi^h would force the processpr to accept T, P, "^Sid 
Q as integer -variables and"<:[^ J, ^ahd \rj as real va^riables for^the.. 
program in whis^h'th^se declarations are^ given. , ' . 



Many otjier languages also eTi5)loy^ forms of ex|)iicit typing. 



Answers ^ to Exerc ises Fg-2 Set B ^ ; . 

2. •'Jg? • and 37 (The space is immaterial.) 



F2-3^ Input - Output Statements ■ • • - '"i^ 

This is tl^e section where t^e student is ,first intrqduced to the concept 

of forjnat. FOFtoAN format seems to b^ a necessary p\±X- No^otlier language 

has. relied a^ heavily on use of format co<Jes as haS FOBTRAN-. Simplified 

input -output 'statements which avoj.d the explicit use of format^ codes have^ 

' unfortunately,' not been available in, mo^t .FORTRA^ lirgplementations^^ T^^ers - ^ 

In the computer fi^ld HaV^, •for this, reason^ long complained of this aspect 

of teaching FORTRAN— but to little avail. On a njore positive note we must*^ 

be (luick*to point out that great fle:?ibility for controlling the precise: forms 

of data input and desired output' is the reward for madjering the use of format 

, f 

code. 



,We teach only a few of its details at first, in order to focus. on the 
main ide^s, ^ • 

Format code is another language- -embedded in the FORTRAN lan^uage--and 
this aifpears to be +be maliX teac^Hrng diffijculty. Format language is us^d to . * 
describe information patterns consisting of group^ of characters called 
fields and groups of , fields called records , which are normally cards for 
input or printed lines for output. It is des(rriptive, rather than active. 
Most FOROgUN ^statements are impeVatlTfes . 36 format code is £nXpplied in a 
declaration which is then .consulted during the execution of- the imperative 
BEAD 'or PRINT statements. . ' , . ' ' ^ ' ' . 

Enlarged pictures of cards' and printed lines' as vi^al -aids might well 
be employed here to goo^ arivaxitage in classroom exposition. 

E- fields ' ' ' V 

It is iielpful for you to know, but not necessar^ily for the stude^it to be 
concerned vith at this point, that an alterriati/ve form of input and output ^ « 
using. "E" notation is" available. An E- field code like E15.'8' might , be used * 
' to input and/or output Avogadro^s number as shown in Figure'^^TF2-l. 

READ 10, AVGAD - ^ ' 

PRINT 10, AVpAD > ^ 

.10 FORMAT (EI5.8) , . \ . 





J 
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Exercises F2«3 get A ' - f 

• Discu^ion - ' ^ . > 9 

\Vhen assigning these exercises you' may wish .to suggest that -a-Il the 

students use the same data values becfeus^ in exercises of Set.C of this 

^s.ection the students may vant to compute resiilts for later -ccarrparison with 

computer output to be developed with the exercij^s in Sectaon^FS-?. One set 

of suggested data values is given* her^i. ' \ , ' ' - ' ' 

^ * < •* • 

For Exercise 1 'use, T = 3.967 j ^ i - ^ - - - ~ ... . 

h = 20, i^= 3, and J = V ^ ^ - / 

a = 3-0,> b = \.-0, c - -2.5, d = 1.5, and X = 2.0 

.2.0 



2^ « 



Anjpwers 



r. 



3 

\\ • " 

5 " 

In Box 1 



u = 12.'5 and v 



r 10.0, ,s =9.0 and, PHI = 1.11977 (a^csine of ^ 



READ and FORMAT 



READ 10, T ' ' 
10 ^X)RM/\T (FI5.8) 



Data Card Pictilre 



2. 
% 



^ TOAD 11, N, I, J 
/ll FORMAT (3 I 15)' 



3. 



a,b,c,a,x 



.READ 12,/A,B,C,D,X 
FORMAT C5 FI5.8) 



If. 



READ X3jj U, V 
13 FORMAT (2FL5.8) 



> T ^ 


. / 3.?<>7 ( 

1 
1 

Ti ■ 

1 

• i ■ . 'V 


■ 

r 

( 

•„ ( 




'} 


/ «?l . '1 . -■> 

J 

n. 1 i ! 1 


i ■ f 


V 

15 • 45 M 75* 
- * 4r 4 4 4 


1 • !' > 

3. vi b j a 

1 ! 


1 

'hi 




/ \2ji\ 2.j . 

1 ' » [ * 
* 1 

1 ' j 

U [ V 1 

1 ^ 
1 1 . 





ERIC 



115 



. I 



|30 



15.' 



A, Y 



r,s^PHI 



REAB 13, ;A, Y 
13 K)BMAT (9 FI5.8) 




10. 



«^ -READ Ih,^ R, S, PHI 
ll FORMAT (3 FI5.8) 



9. 



s 



ARCSIN ro 



J.U977! . 
I 

PHlT 



Answers to exercises F2-3 ' Set B 

In the following answers the^ format ntnnbers which were chosen are 
purely arbitrary. / 

1. PRINT 51, -A,B,J,K,L - 
, 51 FORMAT (2F15.8, 3115) 

or alternatively, though less desirable, 

■/ 51 'FORMAT (FI5.8, F15>8, 115,, 115, 115) ' 

2. ^ PRINT 52, a',B,C,J,K,L,A1,B1,C1 vJ 

52 FORMAT (JFI5.8, 3115) % ' 



•/ 



i 



3* PRINT 53, IKE, JAK, BAKER, CHARLY, D,B 

. 53 FORMAT^ (211^, to5.8) . 



1^ 



Connnent * ,r « ' 

I» Chapter k we wi^l? expand this subject further and consider sprae. 
♦.more complicated problems. For example, suppose in Exercise 3, the Tist 
cAitained a/ seventh i;tem which is real, like T, 

Nov the format (2I15,>'4f15.8) cannot be reused to print the value'of 
T on the second line. Nor may we writ^ (2115, 5F15r.8),-_ because this 
format yould suggest 7 items per line, each I5 columns wide. This 
corresponds to a line width of I05 spaces. We have, however, agreed to 
limit lipe widths 'to 90 ^,columns iii this exercise, - 
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11^ 



•t - 



, Th^ enSHjH'ecord symb6I, wHidh a slksh will be introduced >iji 

. Chapter Its 'iSer^rmits us to write a. format to describe two or more 
prihted lines differing in format, or two or more input, data* cards differing 
in format. . / ' . i , 

Thus, we will be<^able to handle the proposed extension by^^iting 

^ ^ ' ^ PRINT-* 5^, IKE,- JAK, BAKER, CHARLY, D, B, T 

^ 5^ FORMAT (2115, i^m5,8/rL5.8) 



.' Ai\sveys to exercises F2-3 ' Set C 



In Box ^ 



Z . 



PRINT and gORMAT 



PRINT 20, Z 
20 KiRMAT (FI5.8) 




Appearance of Printed Reeult 
( Actual value ) 



6.1^6700000 



2. 



PRINT 21, L 
21 FORMAT (II5I 



.3. 



PRINT 20, Z' 

20 formAt (F15.8) 



36.50000000 



^ PRINT 20, Q t 
^0 FORMAT (F15,8) 



i+.OOOOOQOO 



' ' ^PRIN**20, X 
^0 FORMAT. (Fi5. 8) 



o.i^9ii^fi85 



6, 



AREA . 



. PRINT 20^" .AREA 
20 FORMAT (-Fl^.S) 



5.&r255xxx 



Accuracy to 8 decimal place's is not worth requiring of the student. ^. 
3 or li- place accuracy on all results is sufficient. . The important thing is 
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12 



that the student realize thal^^ the 515.8 field'code will'cause the Lmplter 
^ to 'print 8 digits of^the result to the right of the decimal- point. ( In \he' 
"neict ^et of exercises we wiir ask the student to verify these hand-computed 

values on the computer. ' ! 

Answers to exercises/ Set A 



1. . To express A^^^ ^ 
(a> ^ABSF(A,^.l-5) 



(b) (a 3^) 0,5 

'(c) /SQRTR('a 3) ' 

(d) ABSF(SQRTF(A ** 3)) 

(e) (AJt±l-5) 



(f) .ABSF{A) **^l'.5^ 
•(g) SQRTF(ABSF(A 3)) 



AESF function is unnecessary becaAsW -•A 
is known to be positive so A^^^ mu^ 
* also be^ positive, 

'0,IC,, but expensive ^computationally, 
rerfUire log-antilog procedure)^ 

This is the best way, (Requires the loast 
amoxmt of computatioti, ) ' 

ABSF function is ^unnecessary, 

O.K., but expensive (requirieS the lo^; \ 
antilog process), . ' , i 

ABSF function unnecessary. * -Also* it; 
requires log-antilog process • ! 

ABSF function unhecessary-^herwisei, 
as good as c. i'- " [ 



l3/2 



2. Only f or'^g. would be s&t'i'sf act ory^'ife express [a|-' . £ is 

preferred because it avoids the log-ant ilog process. Another alt^r 
native would be SQRTF( ABSF( A ) ^3). 



J : 

If 



Answers jo exerciges Fg-^ Set B 



In 



4 ♦ 



;he three incorrect statements b^low we see that the first and third 
can be corrected without ambiguity. The second statement may be corrected 
in one <k two ways, but Neither can be said to be preferable. 

1* ' T = 3* - A should be T = B^(-A) • * 



5. 



,F r-^C/-5 + h may b.e 'corrected as F = C/(~3) + ^ 

or as F ='c/(-3 + Ji) 

Saimilaray 

G(= A + B^CC"^ - f/d), may be corrected as 

G = A+B^(C*f-F)/D) 
or 'G = A+B^(C*(-F/d)) 



( These in general' 
jwill yield differ- 
l exit, results^ > ' 



f These are compu- 
tationally 
equivalent . 



^sveirs to exercises F2-6 

1. i(c) and (d) are invali"^ (mixed mode), 
(c) should read Y = L0GF(SINF(F) ) +^ 22.. 
"~dt'^s?i$urd~reyd = J + 1 



2* I = EXPF(Z ^ 2^(a1+*Z + A3)) ^' 

Avoids log ant i- log for' powers of Z. . 

2 < ^ ► ' ^ 

f Factors out Z from subexpression to save 2 mtiltipli cations. 

3. IPART = V ' 

h. IPART = V ! 

FTEMP = IPART ' " ' - -» 

FPART = V - FTEMP 



5. RINTV =.IOTV . _ , . , 

6. 'They are not the same. The box that's n'eeded Is 



J <-">prdnk(i/k) 
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Answers to exercises F2-7 - 



3. 



•* 

1. ' Label 


statement os: declaration 


X 


REAET 10, ^ 


10 


FORMAT '(F15.8) ^ ; 




Z = 2.5 +T 




i.PRIjST 10, Z . 




1 GO TO 1 ^ 




1" 



LaVel 



1 
-11 



21 



.Statement or declaration 



BEAD 11, N, I, J 
FORMAT (3115) 
L = Wil-l) + J'. 
PRINT 21, L 
FORMAT' ( 115 ) 
Gp TO 1 
■END 



^ote: the PRINT stateme;.t may opti^^ly refer to format number 11 ' 
.Which Kould -be equally suitable hfere, thereby reducing the number of 
fonnats needed in this progra^v. ' If f^rn;^ n vere used Tor ^ing, 
executi6n will tenninate ageer printing because the lisfis'then ^ 

exhausted of items in spite ofthS'fact that the governing' format voiid 
have allowed up to thre,e items ^to be printed per line. The ^me ^'client 
api,lieo to the answers to the remaining programs in th;Ls set. ' 



Labet 


statement or declaration 






1 


W A, B, C, 1 X . 






' '''12 


FORMAT (5^15*8) 




> 




Z. ^ {{A^x +'b) *X + C)*X + D 






PRIOT 20^ Z 






20 ' 


; FORMAT (FI5.8) 
GO TO 1 




f 




END 

r 


• 


• 



1^0 



statement'' or declgration 



READ 13, V 
FORMAT '(2F15--8) 
■Q =-SQRTF((U - i^'.5)*V) . 
PRINT ao, Q 

format (fi5.8) 
'go to 1 

END t 



Statement or dedlaration 



REAP 13, (A, ;Y 
FORMAT (2F15.8) 
-X-^ 2./(Y + A/Y) 
PRIHT 20, X ^ 
FORMAT (FI5.8)' 
GO TO 1^ 
END 



0 0 



*^ Statement or declaration 



READ Ik, R, Sj PHI 
FORMAT (3F15.8) *' 
RSQ = R^R 

AREA^ 3. 1^159/2 • * 

-(S.* SQRTF(RS(^ S*,S) + RsVpHI) * 

PRINT 20, 'AREA' 

FORjMg^^tfl^) 

GO TO 1 * 

END 



The dl^it Ji^ in this* colxamn is called a "continuation cod6/^' ' When a 

.statemerlt or declaration is tob lo;ng to fit conveniently on one lin^ it. 

may be continued on succeeding lines. Each of the succeeding lines, if , 

there are more than one must be marked with a non-zero Sigit.' The _ 
- *^ 

maximuto number . of lines permitted for one statement or declaration is 
ngt quite standard. Most FORTRAN implementation^^ permit ug to 10. You 
can double chepk this limit with a reference manttal for your situation, 
although it's not likely to be an important matter. 



1 



TF2 



7^ NUM2a = pki;ce/2ooo 



8. ^ Label 



c 

1 

• 51 

c 

c 



statement . or declara-j^ien 



CARNIVAL WHEEL MODEL ' 

(LET NS ='THE INTEGER S.) ' 

(let np = the integer ' P. ) 

REAt) 51, NS,, M . ; „ , 
FORMAT (2115) 

NS = (M + NS) - ((M + NSV32V32 
■ "tJCMOTE COLOR POSITION, 
*IrE»', (M + NS) HODIH/G k 
K = (M + NS) - ((M,+ llS)/k)^k ' 
COMPUTE PAYOFF, NP. 
* NP = 20*K - 30 ' 
PRINT 51, NP' ' ' /^-^ 

END ' ' * - 



y 



) : 



* . 



*ERIC 
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*^ Chapter TF3* 

Sl^- " • BRANCHING AND SUBSCRIPTED VARIABLES 

I ' Introduction • . - ► ' • 

^ * — ^- 

. ^ • -^v^en the student has mast^ered the flow chart text for Cha4)teii 3 anH this 
con^)anion.5pRTRAN, he is ful'ly equipped to solve computational^ problei^of " 
' ^^^^"^^^^^y complexity. All •the.basic. programming tools are at his ^ 
^ disposals Much of" the material in subsequent chapters falls in the category 
of important -refirfements, Shorthand techniques, example applications and 
' . spec iar concepts,, Uiat why it is very tempting, upon' completion 'of this ' 
chapter, to tarry and solve a large number of problems. We avoid this temp- 
tation as much as ^possible,, because the refinements to be intr<yduced''in each 
o'suc'c^eding chapter make the solution ofyproblems increasitigly and inter- 

esting.^ We expose 'the student to relatively few applications her^-^^rcises 
involving algoritlims emphasize^analysi^ mainly rather .than synthesis, j The 
shift to synthesis is a gradual process which is accelerated ip Chapter*^. 

" " - .J. :* fundamental ideas of Chapter 3^|-branching and subscripted varia- 

, ^les are- mirrored tether easily^ in -the Chapter. requiring few pew ideas. 

For brknching the IF statement is used in FORTRAN. Although its' form 
. permits three-way^ branching, it is used primarily for two-way decision maUng/ 
:^ The expression being tested is always algebraic rather than relational or 
\i=>^ogical^in^ Consequently the s^adent/" must learn how to" recast-delation * - 

like, ; / • . ' ' ^ . 



■ ^s^: A'< B,, .true or false 

into' a form likfe 



A t B, less than zero or not less than zero ^ 
Graded exei^cisds -are used to give the student the practice he needs, 

Th^ new iW'associated with subacr^Lpteci variables is^ that the programmer 
h|s the responsibility^ anii must declare how mu ch>, memory ^s pace is to be - 
all9qated for. each 'vector or array employed in a FORTRAN program. This is 
done through the DIMENSION declaration'. ^ 
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Outline of Chapter • ^ ^ ' , - . 

F3-1 %e conditional IF statement ,is^ introduced. ' Examples , and exercises 
de'^op the relationship beWeen flow cjiart boxes .describing the ^ 
dee-ision, or branching' action and the^ORTRAN equivalent, 

" Socialise th^ IF statement is in the forin of a three -'way branching 
L tool, arid be9)a\!£e wa-vant first to learn how to program sjlmple 

"(two-way) condition-box6a^.the first task i6 that of 1-^a^ing how 
the IJ statement" may- he used in "two-way decision making. 

In the IF statement, expressions are .tested not for true or false,* 

but for -less than zero>,^equal to z^ro or greater than zero. So the 

student is^aught here how to convert or transliterai/e logical ^ex- 

pressi©ns (true or false-) into -arithmetic .expressions whose va^lues 

may be < 0, = 0, or > 0, * * 

* - • . ' ' ' t ^ 

At the end of this section another topic in FORMAT is* treated,^ ^ 
* - . ' 

The H-field code is introduced to show how 'a fldw chartf 'box' like 



if- 



can be coded \ in FORTRM. H- fields mixed with 'I- or F-^elds in 
format code p^Tuit the output of numerical results to W .interspersed , 
with annotation, comments, 'messages, names, etc. Printer^ car3>iage 
y control is ^hown "Uo be easily Accomplished with one-space^ H-flelds 
like IHQ, tor IHl, ^ * ' S ' 

F3-2 Since the flow chart tex"^ covers the subject of auxiliary variables in 
,a way that carries over to FORTRAN In a straightforward fashion,* this 
section consists cnly of one example program plus exercises. ^ • 

F3-3 The transliteratio^of compound condition "box^s "froQi-'flow chart to 

FORTRAN IF statements *is explained and illustrated, '/,Marly-way condition 

boxes can also be transliterated as & series of IF statements. Occa- 

« 

sionally the t^hree-wajf det:ision power of the IF can be taken advantage^ 
of in such transliterations. 



F3-^ SiiKfe no relational symbols (operators) like * ^ 

are available in FORTRAN II, there is no^need to cons idSr* the question 
of precedence levels for such operations, • • 

The FORTRAN forms for writing subscripts and singly-subscripted 
; variables arejjitroi^ucedv The DIMENSION declaration for allocating 

space for one-and two-dimensionai ai-rays is* explaihed and illustrated* 
The notation for .input-output of vector segments is introduced (with- 
^ out explaining the notation as an implied DO- loop). 

F3-6 Everything said about singly-subscripted variables in F3-5 is .^.epeat^d 
Afor doubly-subs^cr;lpted variables. The input-output array segment 
notation is also introduced (but/ again, 'not explained in terms of 
DO- loop concepts). ' , ^ * «. 

• ,. , • • \ \'' . ' ■ ' 

The student is asked ,to work numerous exercises throughout the chapter^ 
which, involve construction c/f FORTRAl^ programs from corresponding flow charts 
presented as examples or developed as exercises in the flow chart text. 



Anfewers to Exercises F3-I , Set A 



i; 


Statement 20' 


• 




Statement 30 


• 


3. 


Statement ^20 


• 


k.} 


_ • 




5. 






D . 


IF^a; iU, iU, y ^ 






K - > 0 i 






t > 0 ' 




9. 






10. 


A - B <^ 0 ^ 


* • 




X - Y - 0^ 






X - -Y / 0 





Answers to Exercises Fj-1 Set 

1. IF(V) 8, 8 • ^ 
Tr(j) 10, 6, 6 
.IF(K - iO'^O, 3t). 
^IF(X - 9.7) 37, 37, 20 

5. IF(I - J) 16, 15., l6 

6, IF(A + B - C) 5, 5, 9 



Answers to Exercises Fj-1 Set C 
1. 



20 



22 



23 



PRim 20, X, Y 

fORr-lAT(5HaXD=o,Flp.3, ?Hayo=oy 91.0.3)* 
i 

PRINT 21, X ' 

FORMAT( IHO, * FIO , 3 , iSHolSoTIIEdVALUEoOFaX ) 
PRIOT 22,, I, J, X 

F0RMAT(lH0,2l6,5HdXo-o^ FIO.3) ' \ 
PRIHT 23 ' 

F0RMAT{^lHlF0imoSC0REoANDOSE\n3^d)rEARSoAG0) 



Answers to Exercises E3--1 Set D 



1, ' Label ' 



20 



21 
. 5 



Statenen^ 



or aecIarlTT5«' 



READ 20>8><r>0^^)j^ 
FORMAT (4Fl5i8) 
PRI>IT 20>B>C>0>X 
IF (B-^C) 5>5>4 . 
PRI'NT ?I>D 
FORMAT (F20.8) 
STOP 

PRINT "21iX 

STOP 

END* 





READ 20>B>C>DfX 


20 


FOR/^AT (4F15.8) 
PRINT 20>B«C^D,X 
IF (0-C) 4>5>5 


4 


T*=C»8+D*X 
GO TO 6 




T= D-C 


PRI^T 21>T. - 


21 


FORMAT (F2'J*8) 

STOP 

END 



3(a); Laijei * Statement or declaration 



V 





r- 




READ 20f8»C»D»X 


* 


• \ 




PRINT 20»B»C»DfX 






•2-0 


FORMAT <»4F15«8) 








T=B+'C 








,UaB»C*D 








IF {T»T+X»X-Uh 7»7»5 






5 


W=T+U 








PRINT -2i»W 


* 






STOP • ' . 


r 




-7 










y>=X*»8 ' 








'PRI"/^T 2L»W»Y* 








STOP ^ 






2-1 


FORMAT (2F20»8r\ : 


✓ 






END 



27 



TE3-1 




* 

^fb)- Label 


Statement or declaration " » 






^ ; 1 — ■ 






PR INT ?Q f B »C • 0 






FORr'AT j4F15»3) 






I F ( OB+C ) **2.+ X^^X-B^C*0 ) 7 »'7 f ^ 






w=*r+R + p*r*n * ^ . 






PR I NT ? 1 • W ■* ■ 






STOP ' , ' 


' 7 










Y=(B+^)**8 






PRINT 21>W»> -^i^i 






STOP • ^ r' • 


^ 21 




FORMAT f2F?C»8) 


4* 




END . 



1|. Label 


Statement or declaration 




. ^ • 1 


READ 15»JiM,M,. 




15 


FORMAT (3P.5»F20,8') 






IF (M-N> 4>4f3^ 




3 


SUM = J+M 






GO TO 5 - 




t\ 


SUM= J+N 




' 5 


PRLNT 15»J»M»N»SUM 






GO TO 1 






END 


* » * • ' 



Label 



S"t;^tement or declaration 



7 




1 

, 25 

5 

. 26 

7 

27 

6 
2ft 

I- 



PRINT 25iB»C ' ' 

FCfRMAT (2F15*8) ' • * * ' , 

IF* (b') 5>J»5 ./ • 

X = -C/P . " 
PRINT 26»X , 

FORMAT (24H. TH6 ROOT OF B*X+C=C -IS F15#8) 
66 TO 1 / . • .- 

^F ( C ) 7 • 6 » 7 ' ^ <vt 

PRINT 27 ^ • , ' ' . 

.STATEMENTS 6 AN© 7 HAVE f-:MPT.Y OUTPgT, LISTS 
FORf'AT (20H B*X+C = 0 HAS NO ROOT) 
GO 10 1 ' - ; . ' 

PRINT 28 • • / . . , ^ 

FORMAT (36H EVERY REAL ^UMBER , SATISF lES B*X+C 
GO TO 1 ^ ' ^ 

END . . • . ' - ' . • 



= 0) 



^2|f ^ 



•^8 



La'bel 6tatement' or declaration 
"~ I I SUMALL k ,0» * , ^ 
NCOUNT 4= 1 
Z READ 20, T ^ ^ 

- 20 F0RMAT(F10.5) 

" ' SUMALIi"= SUMALL + T 

NCOUNT f NCOUNT + 1 
IP(NgOUKT " 100) 2, 2, 6 
6 PRINT 21, SUMALL 
^21 FORf^AT (lOHc3SUMALIo?o, F12,50 
^?TOP ' ^ ^ 

' * END ' ' 



Label 


^ Statement or declaration 






SUMCUB = 0 . 






NCOUNT =^1 ' 


2 




READ 20, T 


20 




F0RMAT(F10.5) 






SUMCUB = SUMCUB + T^3 






NCOUNT = NCOUNT + 1 






IF(NCOUNT - 100) 2, 2, 6 


6 




PRINT 21, SUMCUB 


.21 




FORMAT( 1C«dSUMCUBQ= o, pf 6 • 5 ) 






STOP 






END ' - ^ 



^Label 




Statement or declaration 






SUMNEG = 0 ^■ 






NCOUNT. = 1 . - - - ^ > 


2 




READ 20, 'T 
FORMAT (FIO.5) 
IF(T) 5, 5 


?o 




. k 




SUMNEG = SUMNEG + T 


5 




NCOUNT = NCOUNT + 1 
I?(NCOUNT - loo) 2, 2, 7 * 






PRINT 21, SUMNEG ' 


21 




FORMAT(lOHoSUMNEto=o, FI2.5) 






STOP - . ' \ , 
END 



TF3-1 



^- Xa"bel Starbemeni or decl^ra^tion 



2 
20 



^ 21 



SUMALL" ='0 
SUMCUB = 0 ' ' 
■SUMNEG = 0 
NCOUNT = i 
READ 20, T 
F0RMAT(F10.5) 
SUMALL = SUMALL 
SUMCUB SUMCUB 

IF(T) 5, 6; 6 

SUMNEG = SUMNEG 
NCOUNT = NCOUNT 



T 

T^-^3 

T ' 
1 



IF(NCOUN* - 100) 2, 2, 8 

PRINT 21, SUMALL, SUMCUB, SUMNEG 

FORMAT(lOHoSUMALLa=o, F12.5, 

10HoSUMCUBo=o, FI6.5. 

10HQSUMNEGo=o, F12.5) 

STOP 
END 



5. Label Statement or declaration 







CUMSUM = 0 






NCOUNT = 1 


2 




READ 20, T 
,F0RMAT(F10.5) 
CUMSUM = CUMSUM + T 


20 








PRINT 21, CUMSUM 

F0RMAT( 10HOCUMULATIVEoSUMo=o^ 

NCOUNT = NCOUNT + 1 

IF (NCOUNT - 100) 2, 2, 7 


21 








STOP 

ENlf " . — - * 
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6. Label Statement or Declaration 



2 
3 

5 
7 

11 



6 
• 9 

10 

12 
13 

16 

101 
102 
103 

104 



BAMTMON OR 'volleyball- 
,1 = 0 

ISTATE = 0/ ~ 
JSCA ^ 0 ' ^ 

JS,CB,= 0 

IF(I - lOO) 3, 12, 3 
READ 101, T 
IF(^) 6, 5, 5 
if(istate) 8, 7, 8, 

ISCA = ISCA + 1 

I = I 4- 1= > 

GO TO 2 > ' 

ISTATE = 0 - > 

GO TO 11 

iEF( ISTATE) 10^ 9, 10 

ISTAfer= 1 • , * ' , 

GO TO 11 

iSCB = ISCB + 1 1 ' 
GO TO 11 . . 

IP(ISCA - ISCB) 14-5^ 13, 14 ( 
PRINT 102, SCA " . . / 
STOP \ . ^ 

'IF(ISCA - ISCB) 16, 16, 15 ' 
PRIUT 103, SCA, SCB 
STOP 

PRINT 104, SCB; SCA - ' - 

STOP ^ ^ ■ ' " 

F0RMAT(F12.5) 

FORMAT( iaHDnEoGAMEc3;*l6 . 4HaALL) 
FORMAT( 14HPLAYERqAoWIKSo, I6, 4HaT0o, l6) 

FORMA!PCl4HPLAlffiR6BDWINSe,ol6,4HaTOo, l6) 
END 



Answers to Exercises F3-2 Set A 
1» Label ' Statement or declaration 



2 

\3 



, 6 

101 



FIBONACCI SEQUENCE USED PRODUCE RANDOM^ NUMBERS 
LTERM = 1 * 
NLT = -0 

I =^1 • • ' " 

IF(I - 111^)^3, 6, 6 * . / 
^irCpPY = LTERM. ' 

LOPPING ^OFF ALL BUT RIGHIM)ST IMREE DIGITS , 
LTE3^'= LTERM + NLT - 1000 X ((LTERM +^n£t)/iOOO) 
NLa\= ICOPY^ • < ^ 1 

I = I.+ i ' - - w 

IF(I^- 17) 3, 4^4 • y ^^:r' 
PRINT-^lOl, I, LTERM / ' ""f, ^ 

GO 10 2 
STOP 

F0RMAT(2l6)^ • ^ 

END ' ^ 



2731 



2. Label 



8 
100 
101 



.statement or declaration 



THE OTOSUM PROBLEM 
.EEAD_100/.JDOLD i ^- 

READ 100, TN?W 
WOSUM = TNEW + TOLD 
PRifiT 104r,.TW0SUM 
TOLD = TNEW 
'I fc: I + 1 

IF(I -100) 3, 3, 8 
STOP 

FORMAT(F15*8) 

IPRMAT( 10HoTW0SUMo= F15 . 8 ) 
END ' • 



1. 



3« Label ^ Statement or decla:ratiqp 



ERLC 



9 

100 
101 



8 
9 

10 

11 

* 13 
100 
101 
102 



THE ^LTSUM PROBLEM , 

READ 100, TOLDER 

READ 100, TOLD ' 

1 = 3 ' 

READ 100, TNEW 

ALTSUM = TOLDER + TNEW " 

TRINT 101^ ALTSUM 

TOLDER = TOLD j 

TOLD =-TNEW ^ ' . 

1 = 1+1., - ' 

IF(I - 100) k, k, 9 

STOP ' . ^ ^ 

FORMAT(F15.8) 

F0RMAT(10HqALTSUMo=o,F15 .8) 
END 



Label Statement or declaration 



'MOVING AVERAGE.. 

SEAp 100, K 

BMD 101, TOLDER 

READ 101, TOLD ^ * 

1 = 3,' - ^ 

READ 101, TNEW. ' 

IF(I -'K) 7, 8, 8 

TOLDER = TOLD 

TOLD = TNEW 

1 = 1+1 . ' ' 

GO TO 5 ■ " 

IF(TNEW - TOLD) 9, 10, 10 
AVERGE = 0*5 * (told + TOLDER) 
GO 'TO 11 ^ 
AVERGE = (TOLDER + TOLD + TNEW)/3. 
PRINT 1025 AVERGE ' 
IF(I - 100) 7, 13, 13 
STOP '.a' 
FORMATC16) . 

FORMAT(F15^8) ^ 
FORMAT( 12HQAVERAGE oES □, FI5 . 8 ) 

END > 
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5/ La-bel 
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3 



6 

101 
102 



Statement or declaration 



THE REGIONS TABLE* 

PRINT -101-^- — ' " 

N = O ' ' 

lA = 1 

ISUMA = 1 

IB = 1 . 

ISUMB = 1 

IC = i 

PRINT 102^ N, lA, -IB^- IC 

- 15) 5, 5, 6 # ' 
IC = ISUMB + 1 
ilB = ISUMA + 1 
^lA = lA + L 
|J = N + 1 ' ' ' 

SUMB = ISUMB + IB 

SUMA = I^UMA + lA 
CO TO 3 

stop 

FORMAT( 2l^HDG0DD ifooQo qAdddod BDoaooC ) 
F(DRMAT(4l6) 



) 



Answer to,. Exercise F3.2 Set B 

-Label Stattement or declaration 



5 
6 
11 



7 
8 

• 9 

10 

100 
101 



REAB 100, KG, KD 
KA 4 KG 
ICB 4 KD 

IP(KA - KB) 5, 5, h 
KR = KA 

GcTtQ 11 

IF(KaI) 6, 7, 6 

KR = Kb - kb/ka*ka ^ ' 

KB == KA 
KA t:> 
GO TO 

IP(KB)i9, 8, 9 
KX = 0] 
GO TO m 

= K*KD/KB 
PRINT ibl, KC, KD, KX 
S!SOP ' ! 
F0RMAT(^I6) 

formatJ in 0,316) 

END' 
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— Ans^fers to Exercises F3-2 Set C 



1, Label 



Statement or declaration 



1 

20 
21 



READ 20»X1»Y1,X2»Y,2 
PRINT 20»Xl»YlfX^»Y2 
FORt^AT (4F15,»3) 

DI5T= SQRTF( ( X2-X U **2+ ( YJ-Y 1 ) *.*2 ) 
PRINT 21/DIST 

FORMAT (21H THF LENGTH OF PQ IS Fl5< 
GO Ta 1 * 
liE'ND ' 



3) 



2, Label 



Statement or declara-edon 



1 READ 20'.Xl»Yl,X2fY2 
'PRINT ^r),XlfYl,X2»Y2 * 
20. FORMAT (4F15»3)* 
^^'*T IF (X2-X1) 4t6»4 
"4 S=(Y2-Y11/^X2-XH 

PRINT 21»S , • , ' • , . 

21 FORMAT (20^^ THE SLOPE OF PQ IS FfiOtS) 
GO TO 1 ' 

22 FORMAT (29H PQ IS PARALLEL TO THE Y-AXIS) 

GO TO 1 ; 

lENb • ^ 



3. Label' Statement or declaration 



1' 
20 



21 



8 
9 
22 



10 
23 



READ" 20»X1»Y1#X'2kY2 

PRINT 20fXl»Yl,X2»Y?^ 

FORMAT (4F15#3) . ' U 

READ 20»DELX * 

PRINT 2h, DELX 

FORMAT(8HoDEIXD=c(F15.3) ' ' ' ' . J 

IF fX2-Xl) 5»8,5 - 

S= (Y2-Yl)/(X2-i(l) . , 

DELY5S*DELX I 

PRINT LltDELY 

FORNAT (8H DELY = F15»3) 

GO TO 1 

IF^ (DELX) 9»10,9 
PRtNT 22, ' , 

FORMAT t21H NO SUCH VALUE EXfSTS) 
GO TO 1 

PRINT 23 ^ > ^ ' 4 

FORMAT (24H AMY REAL NUMBER WILL DO)-^ 
GO TO 1 ' 
END 



ERIC 




3034 



Label Statement or declaration 



1 
20 

% 2k 
* 5 

8 
21 



IP 
12 

.23 

11 
22 



READ 20fXlfYlfX2»Y2 
PRINT 20f X1»Y1,X2»Y2 
FORMAT <4F15«3) 
REAr 2C»DELY 
VRim 2k, DELY 
EORliAT( 8HoDELYo= qF15 . 3 ) 
IF (Y1-Y2) '5,10,5 
IF (xr-^2) 6',9,6 
S= (Y2-Y1 )/(X2-Xl) 
D;ELX« DELY/S " 
PRINT 21,DELX 
FORMAT (8H DELX = F15^ 
GO TO 1 ' • 

DELX=0 
GO TO 8 

IF (DELY) 11,12,11 ' . 

PRINT 23 

FORMAT (24H-'ANy REAL • NUMRpR WILL DO) 
GO t;o 1. ^ , ' 

PRINT 22 

FORMAT (21!I NO SUCH VAIUE EXISTS) 
GO TO 1 

END. • r , 



»3) 



^> Label Statement or declaration 



20 

2k 
-^5 

21* 



8 
9 
22 



23 



^1 



T 



READ 20,X1 ♦Y1,X2,Y2 
PRINT 20*,x:,Yl,X2,Y'e 
F.ORMA^T (4Fia»3-) - ^ 
READ 20^X? _ - , 
PRIIJT 2k, X ' r t- ^ 
FORMAT ( 5HD353=:o>l5 . 3 ) 
IF rxi-X2) 5^8', 5 
S= 'tYZ-Yl) /(X2-Xi:) 
f= YL+S*(X-xr) ^ 
PRINT 21, Y 

FORMAT (5h 'Y^=. F45.3) 
30 TO 1 

IF (X-Xl) 9,10,9 \ 4 

^RTNT^TZ""^ 

FORMAT (21H NO SUCH VALUE EXISTS) 
GO TO 1 * " .'^ 

PRINT* 23 y , 

FORMAT (24H ANy' REAL NUMBER WILL DO) 

60 TO 1 

END 



er|c ■ _ . . 



-J 



^ Statement or declaration 



1 REAP"?OfXl'»Yl»X2»Y2 . ^ 

PR>INT 20fXl*Yl»X2»Y2 ^ 
20 -FORMAT (AFTSt?) 
READ 20*Y ' 

PRINT 2lt, Y • ' 
2lf FORMAT(5HQYa=rQF15,3) 
IF;(Y1-Y2) 5fl0#5 
5 IF (Xi-X^) 6f9f6f 
•6 S=(Y2-Y1) A(X2-X1) 

X= X1+ <Y-Y1)/S 
' 8-PL£LWrr 21»X 

21 FORMAT (5H X = F15»3) J 
60 TO 1 ' . 

9 k~Xl 
• 00 TO' 8 

10 JF (Y^Yl) llfl2fll ^ 

11 PRINT 22 ' • ' ^ 
22,FORhAT (?1H NO SUCH VALUE ' £X I STS T 

, 60 TO 1- ^ 

12'PRINT 23 \ . , 

23 FORMAT' (2^H ANY REAL NUMBER WILL DCT) 

60 TO. 1 , . 

END 



mil 



' Label 



Statement or declaration 



•20 
9 



/ 

11 



22 
*23 
24 



RE/fD 20>XI fYlfXZfYZ 
PRINX •20>XljYl»X2»Y2 
FORMAT {4F15»3) 
IF (Xi-X2) 3>lli'3 
IF (Y1-Y2)' 4f9»4 
PRINT 21 
PRINT 22»Y1 
60 TO .1 

S:^(Y2-Yl)/rX2-)^l) 
XINTa Xl-Yl/S 
YINT.= -S*XOT 
PR-INT 23»XINT " 
PRINT 22fYfNT 
60 TO 1 ' 
PRINT 23fXr 
PRINT 24 

60 T6 1 _ 
£0RMX^-~(-29H "Pd DOES •NOT 
FORMAT (16H ^-INTERCEPT 
FORf^tAT '(16H X-INTERCEPT 
FORMAJ (29H Pp DOES NOT 
ENq 



INTERSECT X-AXIS) 
IS F15#3) . 
IS F15.3) 
INTERSECT Y-AXIS) 

> 
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: 8, Label 



Statement or declaration 



20 
I 
3 



6 

10 



7 
8 

9 
12 
11 

ai 

22 

23 
24 



V 



o 



^READ 2'6,X1,Y1,X2,Y2 
PRINT 20.X1,Y1,X2,Y2' 
FORMAT (4FI5.3) 
IF (XI 3,11,3 
S = (Y2-Y1)/(X2>X1) * 

nm = xi-Yi/s 

YIKT = Y1-S*X1 
IF (Y1*Y2) 7,10,10 
PRIOT 21 
GO TO -8 
PRINT 23, XIKT 
IF (X1*X2) 9,ia,12 
PRINT 22, YINT 
GO TO^l " 

PRINJ 24 ' - — - ^ 
GO TO 1 ^ 

XINT = -^Xl' " - 

GO TO 6 ^ 

FORMAT ^9H PQ DOES 'not INTERSECT X-AXIS) 
FORMAT (I6H Y-INTERCEPT IS-FIJ^) 
FORMAT (i6H X-IKTERCEPT IS FI5.3) 
FORMAT (29H PQ IX)ES NOIJ^JNTERSECT *Y-AXIS) 
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33 
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Answers to Exercises F3-3 Set A 



1« Label Statement or declaration 



IF 



-(2.0-X) 
(X-7t0) 



2f 2i30 
i'20f 20f.30 



^2. Label 



Statement or declaration 



IF (7.C-Q) 20f2»2 
IF (7#C-R) ^f3f3 
IF (7»0~S) 20f30»30 



3« L^bel Statement or declaration 



IF (1#7-X) 2»30f30 

IF (X-8/4') 3f30f30 

IF (-3#9-Yl 4»30f30 

IF (Y-5#4) 2'0»30>30 



^Ca) La^el Statement ot d^laration- 



IF 
fF 
IF 



(XI ) 3C.»30f2 
(0#5»X1-Y1 L 3f 30*30 
(Y1^-2^*X1) 20»30»30 



k(h) Label Statement orj^eclaration 
ifiOf 3 



IF (XI) 5;30T30 
IF (2#0»X1-Y1) 6^30»30 
IF (YW'^tS^XlJ 20»3'0'V30 



kic) Label Statemenj: or deciaration 



IF (XI) 5,'30»2 r 
IF% (0#5*Xl-ri) 3»30»30 
IF (Yl-2#r"'X.l)'^a»30»30 
IF (2»O»Xl-Yl)-6^'30y30 
IF '(Yl-0.5^^Xi/ 2a»30»^O 



^. Label Statement:; bjr*"declaration 



IF (XI >^:3a»30»2 
IR (Yl ) 3b»3Crf 3 

\f (Yl-t-r.6667*Xl+2#^ ) ) 20t^0f30 



6. Label Statement or declaration 



IF (XI-3.li+l59) 2,2130 

IF (C.5^(3.141^9-X1)-71 )' 3f^f^0 

klF fYl--SINF(Xl -) )r-20f 20f 3Q- . 



1 



7. Label 7" Statement o'biiq^jaration 



IF (Yl) 30,30,2 

IF a(Y1 - {'h.O * XI 

IF (Yl - lh,Q * n - 



H 16.0))30,3,3 
12.0))30,20,20 
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> Answers ^ to Exercises F3-3 ySet B - 

!• ^ ^ IF(I - 7) 2, 10, 10 

10 IF(I - 20) 3, 3/^ 

2. ' IFCW - A) 1, 10, 10 - 

10' . ■lF(W (B - 2)) 2, '11, 11 

11 IF(W - C) 3, 3, h 



3, Label Statement or declaration 



13 

10 
11 

12 
✓ 

20 
21 
22 
23 
24 
25 
26 
27 



READ 20»Xl»Yl»X2fY2 
PRINT 20»X1»Y1»X2»Y2 

IF ixi) m^»7 

IF(Yl) 9»6»8 
PRINT 21 
60' TO 1 
PRINT 26 
-TO 1 ' 
PRINT 22 ^ ^ ^ , 
Gb TO 1 
PRfNT "25' " ; 
IF (Y^) laSfl 
PRINT, 27 
GO. TO 1 

IF (YD 12»6»11 ' 
PRINT 23 " , ^ 
GO Ta' 1 ' , 
PRINT ?4 
GO TO 1 , 
FORMAT (4F15*3) 
FORMAT (2H 1) 
FORMAT (2H 4) 
FORMAT . (2H 2) 
FORMAT (2H 3) 
FORMAT (21H P 
FORMAT (21H P 
FORMAT (35H P 
END / 



LIES ON TrfE Y-^AXIS) 
LIES ON TfiE. X-AXIS.) 
1S THE ORIGIN) 



ERIC 



39. ' 



35. 



Label 



Statement or d'eelaration 



-1 
51 

40 
41 
42 
5 

6 

7 

8 
' 9 



Carnival wheel mo»€l 

READ 51»NS»M ^ 
FORr^AT- (2P5) 
K = M+NS.- ((M+NS)/if)*if 
IF (K) 41t''^»41 
IFU-1) 42»6»42 
IF (<-2) 7»7,8 

GO TO 9 
NP= "-30 
GO Ta 9 * 
NP = 0 
GO TO 9 
NP = 50 

PRINT 31»NP 

GO TO 1 f 



Comment : For more efficiency we^ could use just twcJ IF statements. 
Statement hO w^uld b'e removed. Statements hi and k2 would be 

hi iF(>K - 1) 5, 6, i|2 
1|2 IF(K -'2) 7, 7, 8 ' 

We use the first IF'^as 'a-three-way branch here. 
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F3-^ Precedence Levels for ' Relatio ns ^ ' . 

- — ' , _ • • 

In more advanced versions of FORTRAN, especially FORTRAN TV or variations 
of i^. relational operators can be used in an IF statement called' the "logical 
IF" as dist languished from/the one discussed in this chapter which is usually* 
.called 'the "arithmetic IF." " • 

^ ^In the/fogical IF one could write for example * . 



which is equivalent to 




Here the character group ^kLT. means "is le^s.than," Other character 
groupings are used to^ designate the ^ive'ot^er relational operator symbols.- 
all 'of which are listed here. - ' 





,< 


< 








> < 






.LE. . 






:gt. 





As another example^rrCSne 'could write ' . 

which is eql^Lvalent^to . ■ . ' \- 




NEXT 



ERLC 



Precedence Revels for these relational operators are exactly 4s -described 
in Section- "'^ ^ ^ ^ - 



Ansvers to Exercises F3-3 *Set A • . . * . 

1. ^X(5) * . ' . ' ^ • ' 

2. X(J) ^ * ' , * ^ 

3. CHAR(l) ' use upper case for letter i. ^ , > . 
1*. ' B(I + 2) ^ , ' 

5. X cannot be used as a subscript since' it is not an^pteger variable, 

6. Z(5*I + 2) , . . ' \[ , , 

7. I + y is an invalid torm for a subscript expression f variable 1 +, 
variable 2 not allowed)^ . . . ' * j 

- 8. not valid; only unsigned integer constants (non-zero) are permitted; 

9. xyz(nY) . . 

10. -5^J would be an invalid subscript expression. An unsigned integer 
constant (non-zero) is called for in this-. 



Ansvers to Exercises Set B 

'1. DIMENSION A(50) 

- READ 5, (A(l), I = 1, K) 



' 2. DIMENSION S(l25) ' . 

READ 6^ (B(J), J = 5, -N, 2) f ' 

3. DIMENSION A(50), B(50) 

READ (A(T), I = 1, lO), (B(i), I = 10, N, 2))/ ^' 



1. » Label Statement or declaration 



n 

51 



CARNIVAL WHEEL WITK SUBSCRIPTS 
DIMENSION NP(tf) 
,READ 51»(NP(I)»I=1»4) 
READ 51 »NS»M 
FORMAT (21 IS) 

K=M+NS - ((M+HS)/4)* 4 
PRINT 5UNP(K+1) 
GO TO ip 
END 



2. Label S-^atement or declaration 



>20 



21 



22 



7 

10 
23 
11 



FIGURE 3-29 
DIMENSION B( 100) 
1 = 1 • ' *^ V'f 

IANY=0 
^^AD^pK 
FORMfr(F10»5) 
READ 21»(B( I)»I = 1'»100.) 
FORMAT(7F10.5) 
IF (B(n-C) 5»8»8 
IANY=1 

PRINT ZZfUB^i) 
FORMATC I5»^'10.5) 
I + l 

IF (I-lbO) 4»4»7 ' 
IF (lANY) 11»10»11 
PRINT 23 . > , 

FORMAT (5K NONE) 
STOP 
END 



43 

39 



Label Statement or declaration 



100 
101 
102 
103 

ro4 



FIiroiNG THE ACTUAL DEGREE OF A POL^CMEAL 
DHASION A (51) 
READ IpO, N 
M = N + 1 

m IS FORTRAN UPPER* BOUND (ONE HIGHER) 
READ 101, (A(I), I =.1, NN) / 
ifCa(nn)) 6, k, 6 
NN = NN - 1 

LOWER BOUND IS ONE HIGHER ALSO 

IF (NN - 1)^7, 3, 3 

PRINT 102 

GO TO J. 

N = 'NN - 1 

PRINT 103, N 

PRINT 1(A, (A(I), I = 1, NN) 
GO TO 1 • 
FORMAT (15) 
JORMAT(5F10.5) 

FORMAT( 20HaP0LYN0MIALarSdEMPTY) 
F0RMAT(22HpP0LYN0MIALaDEGREEarSO, 15, 
31H.DDC0EFFICIENTSaAREaASaF0LL0WS. ) 

FORMAT(lHa,5'F15.5) 
END 



Answers to Exercise F3-5 Set D 

Lab^l Statement or declaration 



1 
100 



10 

11 

1€ 



101 



ORCHESTRAVILLE 

DIMENSION IAf5Q0) 

READ 100/ N, (IA(K),K = 1',N) 

F0RI^AT(10I5) ^ - . 

FIRST DATA CARD HAS VALUE OF N IN THE FIRST FIELD, 
THEN UP TO 9 AGES. 
K = 1 . 

IF(IA(K) - IA(K+L)) .6, 6, 5 
K = K-fl 

IF(K - N) 1+, 8", 8 , 
ICOPY = IA(K) 
IA(K) IA(K+1) . 
IA(K+1) = ICOFY 
GO, TO 3 

MID = N/2 ^ . , ^ 1 ' / . 

TEST N TO SEE IF IT IS EVEI^ FOR COMPUTING ! MAD 
IF(MIDi^2 - N) 11, 10, 11 > ' • \ 

MAD = MID + 1 ' : 
GO TO 12' . ^ ^' 
MAD = MID / • i ' 

SUM = IA(MID) + IA(MAD) ' 
TMEED = .5 * SUM' ' ' 

PRINT 101, IMBED, IA(i),IA(N) ' , . 
FORMAT ( l5HcmDIANfclAGEaISOl5 , 
13HaY0UNGESTdrSO, 15 , 11HD0LDESTI3ISO, I5 ) 

GO TO 1 * ^' ; 



Note that several name^ of variables must be changed ^from ^ 
the flow chart in order to 4es-ignate' the desired integer or., f ioa^t^Sg, 



modes* 



ERLC 



ho 



used in^ 
jppj-nt , 



■'44 • 



i ■ - 



1^ • 



Answers to* Exercises F3-6 



1- Lal)er Statement or declaration 



DIMENSION P(22, 2?) 
COLSUM = 0 , 
1 = 1, 

IF(I - 12) k, 5, h 

COLSUM = COLSUM + P(l- K) 

JF(I . 22) 6, 7, 7 

1 = 1 + 1- 

GO TO 3 ^ 

PRIOT 50, COLSUM^ * . . 



Label Statement or declaration 



DIMENSION P(22, 27) 
J = 1 

P(L, J) = P(L, J). + 
IF(J - 27) 5, % 

= J ^ 1 
GO TO 2 



P(M,J)' 



3* Label -Statement or declaration 



DIMENSION P(22, 27)' 
J = 1 

IF(J - K) 3,.-^, 3 

P(L; J) = P(L, J) +2.0 X'P(M, j) 

IF(J - 27) 5, 6, 6 

J = J + 1 

GO TO* 2 
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h, ^ Lat^el Statement or de.claration^ 



DIMENSION ^(22,27) 
J = 1 

COPY = P(L,J} 
P(L,J) = P(M,J) 
P(M,J) = COPY 
IF(J - 27) h, 5, 5 
J = J + 1 ^ 
GO OX) 2 



Label Statement or declaration 



( 



10 



DIMMSION P( 22,27) 
1 = 3. ^ 



MAX = 0 

IF(ABSF(P(L.J)) - ABSF(MAX))5,5,^ 

MAX = P(L.J) 

IF(J - 27) 6, 7, 7 

J =• J + 1 

GO TO 3 

J .= 1 

P(L,J) = P(L,J)/MAX' 
rF(j - 27) lOy 11, 11 . 
J = J + 1 
GO TO 8 



\ 
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Chapter TFk ' 
LOOPING 

* • 

The main sections of this chapter are; 

Fl^-l The DO statement * , ^ | 

Illustrative examples ^ ? 

Table- look-up 
Fk'k Nested DO loops 

^ This chapter' follows cl-osely Chapt.er k of the Flow Chart text. 

- ■ 1 

Fl^rl The DO statement is FORTRAN'S somewhat restricted counterpart, of - 

the iteration box. It is introduced ia this section and many of its 
limitations are spelled out. to the student. The notion of the range 
_ or scope of an iteration or DO- loop is-^tressed and. th'e CONTINUE ■ 
statement is explained for marking the end of the sQope or range of the 
DO loop. . 



You shou^ji check -tha refetence manual for the FORTRAN implementation 
which you will be using and determine if it is one which always causes • 

execution of the scope at least once (most common type) or one whic; 

(like the -iteration box) checks first to see if the DO- loop variable's 
initial value has r^ot already exceeded the "final" value. Most lEM 
processor^ for FOR!I?RAN II fall in the first category. Nearly all CDC ' 
processors for FORTRAN II fall in t^e second category, etc. . » 

^^"^^"'-S^^''^™^^^^ iteration.box used for 'simple^ lopps in 




ectidR-2 are mirrored in thts section, using DO-loops. Much attention 
must be paid ^6 the, problem of getting the student into 'the habit of . ' 
having all statements witjiin a loop ultimately lead to the CONTINUE 
statement of that loop for guaranteeing repetition. The pitfalls are 
described fn c/onrlection with Figures ?k-^ and Fk-C. 

The. last exercise of^^et C brings up a point about the naming of 
library functions which was not touched on when these fiinftions were 
first , introduced in , Table F2-2, ^ - ' \ 

it 

The ABSF is useful fo^r both real and integer arguments. But ABSP ' 



f;or instance is defined for real arguments only. ABSF(x) for instance 
returns a real value. Occasionally we may want to express the absolute 
value of an integer "variable or expression. To do this we can use the 
ABSF function provided the letter X is prefixed to the name/ as 
XABSF. 

Fi+-3 Th^ main purpose of this ^section, beyond giving the student more 

practice in converting 'flow charts having iterations into correspond4nig , 
FORTRAN programs, is to treat certain ' details on input-output. ^ 
These are 

1. The use of implied DO- loop notation as, list elements of input or 
output statements. This notation was used in Chapter 3 without^ 

» explanation there. Here they are explained as a special shorthand 
related to DO statements. 

2. The X- field code for format .is introduce*' and illustrated."^ 

Flf-lf - Details on the nesting of DO loops are discussed by taking some of 
4 the flow charts, in ^ Section \'\ and showing how these may be coded in 
FORTRAN. The doubly -fmp^i^d DO-loop notation for input-output list 
elementa is explained in terms of nested DO-loops. This discussion is 
daggered to suggest its urtimportance relative to the other topics that 
are covered. » » 
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Answers 
t 



to Exercises 



Label Statement or declaration 



60 



70 
50 



DO 50 I = 1, 50 

^feEAD 60, ID, A, B. C 

FORMAT(llp,3F10.5) 

D = SQBTF'U*^ + B^^2 + 

FBim 70, ID, A,B,C,D ' 

FORMAT( IHa, 110 , 1*F10 . 5 ) 

CONTmjE 

STOP 

END 



Label Statement or declaration 



100 
5 

60 



70 
50 

80 



REA£) 5, N 

FORMAT (15) 

DO 50 I = 1, N 

READ 60, ID, A,B,C 

F0RMAT(I10,3F10.5) - 

D = SQRTF(A**2 + + C^2) 

PRINT 70, ID, ■A,B,C,D 

FORMAT( IHa, ip, 4F10 .5 ) ' 

CONTINUE 

HRINT 80 

FORMAT( 13Ht3ENDoOPtiTABLE ) 

GO TO 100 

END 



Label Statement or declaration 



40 



30 



LTERii = 2 
NLT*= l''^ 
IS = 1 

DO 30 I = 1, So 
PRINT 40 I, LTERM, IS 
F6HMAT(lHa,3ll5E) 
IS = NLT + IS 
ICOPy = LTERM , 
LTERM = LTERM' 

ICiT = icopy 

CONTINUE 
STOP 
END 



■\ id- 

ERIC • " ' 



Alternate solu^bo^^wi subscripts, (Note the trick of shifting index 
values to^avpid using a zero subscript ±n FORTRAN,) 



Label .Statement or declaration 



'roi 

60 



bimehsion if(63),is(63) ^ 

IF(2) 1 
IB.= 0 

isdi ^ 0 r 

DO 60 I = 1,62 

IS(I) = IF(I) + IS(I-I) 

IF(I+1) ^ IF(X) + IB \'\ 

IB = IF(I> 

IF(I - 3) 60, 5, 5 

TO PRINT FLOW ^IHART gUBSCRIPT VALUE 

IMl =T I - 1 ' 

PRINT 101, IMl, IF(l)i IS(l-2) 

KiRMAT(lHa,3I10) 

CONTTNUE ^ 

STOP - 

END ' 



LalDel* Statement or declaration 



51' 



20 
.30 



i 



CARNIVAL WHEEL PROBLEM REPEATED 
DIMENSION NP(i^) . . ' 
READ 51, (NP(I), I = 1,4) 
READ 51, N 
FORMAT (^^115) 
ISUM*= 0 
DO 20*L = 1, N 
READ 51, NS, M 

K = M+NS - ((M + NS)A)*J^ M 
ISUM = ISUM + NP(K + ) . 
CONTINUE . ' 

PRINT 30, N, ISUM 
F0RMAT(7H AFTERO,I15, 

30a3PIjyS,aYOUH3NEraWINNINGaaAREa I15,7HdPOINTS) 
STOP ^ „ ' > . , ^ * 

END . ' . ' 



k6 ;- / 

50/ 



5^ --La^el Statement or declaratlnn 



C 
C 



101 



102 



103 



50 
105 



PAYROLL PROBLEM. MAXIMUM NUMBER 

OF EMPLOYEES ASSUMED TO BE LESS THAN 1000 

DIMENSION T( 1000), R(lOOO) 

READ 101, N ^ 

F0RMAT(I5) 
PAYROL = 0. 
READ 102, (T(I), I 
ONE TIME PER* CARD. 
P0RMAT(F10,2) 
READ 103, (R(I), I 



= 1, N) • 

IN HOURS AND HUNDREDTHS 



, , - = 1> N) 

UP TO If RATES PER CARD/ DOLLARS AND CENTS 
PER 'HOUR 
FORMAT (too. 2) 
DO 50 I = 1, N 

WAGES = R(l)*T(l) 
PAYROL = PAYROL + WAGE^ 
PRINT 10k, I, WAGES 
FORMAT(lHq, 110, F15 • 2) 
COIJTINUE 
PRINT IC^, PAYROLL 
FpRMAT(lHD,^25.2) / 
STOP / 
END • / 



7 



51 



Answers to Exercises Set A . ; . 

1* Statement 2 should read 

» ' ' DO 4 J = 2, N 
was missing 

ai^o the statement^following the IF shovdd'read 
) 3 MAX = absf(a(j)) 

label was missing * ^ , . 

'2. DO statement should read 

PO 10 J = 4, N,2 . ^ 
no comma allowed here! 

r 

Change MAX, wherever it appears, to a real variable like MAX. 

# 

3, Change FACT to an integer variable, like IFACT or NFACT, 
The DO statement is incorrect as written. It should read ^ 

, ' DO 10 K = 1, N, 1 
or DO 10 K = 1, N 



k. Either the DO statement should read 




DO 10 K =r 1,- 

♦ 

not L 

or the PRIOT -statement should read 
PRIOT 5VIi, LTEBM 

n 

not K • • 

,also copy should be replaced with an integer variable like*^ ICOFT. 



'TPll. 



N Comment on the discussion which immediately foUoW Exercises F^^-2 Set A. 
- -Oftin we use an iteration box to control the manipulation of a set of 

. elements where the nMtnber of elements, represented by N, is previously 
assigned ja value. Sometimes we would like to assign a value to ' N to suggest 
the en^tjl set so that when the iteration box is encountered the initial test 
is false, with the result that traniiKs through the 16op are taken (i.e., 
no set elements are majiipulated)* '' ' 

A good example arises in connection with a problem to be solved "in 

Section 5-lt {Set C Problem 2). -mere we are to write a procedure called' ' . 

SIMPLIFY, to simplify the integer .coefficients of an arbitrary pol^ncsnial 

given its degree n and ■ the coefficients a' a . ... a ' 

0 1' ' .n* 

Siiiiplificatioh amounts to' dividing each coefficient by the absolute value' 
of the largesj; comon factor. However, when the pol;ynomial is identicall;y ' 
zero, and this is implied whe« n has been assigned a negative value like 
-1, our iteration box allows us to automatic all;y avoii^ the disaster of dividing 
by jzero. JChus in the loop^ shown b^w: { ' ' 

L_ 

i < n 



■ 0 -- 



i«-i+l 




^ivide a^ 
by largest 
common factor 



if , n - -1 at the time we enter Box 1, we will g^ immediately.fd Box 2 with- 
out ever entering Box 3,. because 0 < -1 is false. 



^ ^ Nqy, when coding this particular loop in^FORTRAN there are two problems. 
First,- the n + 1 integer coefficients of the polynomial 

Sq' V ^2' ^n , , " ' 

cannot be assigned to variables IA(o), IA(l), IA(2), iA(n), ..because 

zero subscripts arenot allowed i;i FORTRAN. So we must shift the index of ' 
all the coefficients; That is-, • a^ is now associated with, IA(l), a is* 
associated kLth'^(2), etc. Ihis mean! that the initial value for the DO 
loop index I is 1 instead of zero and the upper limit nrnst have the yalue 



n + 1 . instead of -n.' That^s easy to fix. We dan write, say 



NP1,= N + 1 
^ 50 I = 1,NP1 - 



Now we notice 1;he second problem. If N is initially assigned -1, then 
for the first .value, of I = ^/x^^ exceed^ UPl which is now' 0. If we 

are using a FORTRAN II processor which generates code that bypasses this 



test on the first go aroiiMy we now need tq include the type of test mentioned 



in the st^udent text. For example: 



IF(N) 3, 3, 1 . 

NPi = N +' rap 

DO 50 I = l,^Pl 



\ 



V 



ERIC 



50 54 



Answers to Exercises^ Set 6* 

1. Label Statement 







IX> 10 I = 1, N 

COPY = P(I) 




> 






'P(I) = Q(I) . 
Q(l) = COPY 






CONTINUE » * 



2. 



Label Statement 







DO 10 I =^2, N 2 






COPY = P(I) 






P(I) = ^(I) 






Q(I) « COPY 


" 10 




CONTINUE 



I 



Label Statement 



DO 10 I = 5, N, 3 
COPY = P(I) 

P(I) = Q(I)- - 
Qfl) = COPY 
10} CONTINUE 



Label Statement 



DO 10 I = 1, ND2 

Q(r) = P(I) - 



,Labe;i Statement 



10 



N02 •= 1^2 

DO 10 I = 1^ N02 

J = N02 + I \ 

Q(i) = p(j) * 

CONTINUE 



The altemate solution suggested in 
the' Teachers* Commefttaty is, not ; . 
possible in FORTRAN. You cannot write 
DO 10 I = 1, N^ 



This is illegal. 



^ '^'^ 



■Note: You cannot write P(N02+I) in FORTRAN. 



invalid subscript 'form 



6. 



7. 



Label 


Statement 






N02B = N/2 






IF(N - (N/2) ^ 2) 3, 6> 3 


^3 




K = N02B + 1 






GO TO 1+ 


6 




K = 'Np2B 






DO 10:.I,= 1, N02B' 






KK =■ K + I ' ' . 






Q(I) = P(KK) 


10 




comNUE 


Label 


Statement 



10 



DO 10 I 
LLL V N 



= 1> 
+ 3 

L = N + 1 - I 
P(LLL) = p(l) 
COUTINUE 



K 



'FOR^TRAN doesn^t allow us to decrease 
th^'-value of the loop 'eoun1>er as was 
dpne in- the equivalent flow chart box. 



8(a). 



8(b). 



Label 


Statement ^ * , 








SUMCUB = 0. 








DO 10 I = 1, 100 








SUMCUB = srJMCUB + 


P(I)*^ 


10 




COUTINUE 








SUMNEG - 0. 








DO 10 I = 1, 100 








IF(P(I))1+, JO, 10 








SUMNEG = SUMNEG + 


p(i) \r 


10 




CONTINUE ' ' 





8(c), 



5 
6 

7 

10 



SUMCUB =0. ^ 

SUMNEG = 0. ^ 

SUMBIG = 0. 

DO 10 I = 1,^00 

SUMCUB = SUMCUB + P(l)^^ 4 

IF(P(I))5, 6, 6 

SUMNEG = SUMNEG + P(l) 

IF(ABSF(P(I)) - 50-) lO, 10. 7 

SUMBIG = SimiG +^ABgF(P(l)) 

CONTINUE 



9»' Label Statement 







COLSUM = 0. 






jDO'lO^I = 1, 22 






IF(I -12) 10, h 
COLSUM = COLSUM^ P(I,K) 


' It 




10 




CONTINUE 






fRINT 15, COLSUM 



ERIC 



5G 



TFk 



10. Label ^ Statement 



10 



DO *10 J = 1^ 27 

P(L,J) = P(L,J) + P(M,J) 

COmtfNUE 



11, 'Label . Statement 



3 
10 



DO' 10 J = 1, 27 

IF(J - K) 3, 10, 3 

P(L,J^=: P(L,J) + 2,0 ^ P(m;j) 

CONTINUE ^ 



> 12 . Label Stat ement 



10 

h 

3 
5 



DO 10 I = 1 N 

IF(ABSF(P(I)) - 50.) 10, 10, 3 
CONTINUE 



lANY = 0 
GO TO 5 

W = p(i) 

lANY = 1 



This Is a, little more tricky to coSle 
in FORTRAN than "it may appear at, first. 



13 . ^Label Statement 



10 



DO 10 I = 1, N 
J = N - L + 1 

IF(ABSF(P(J)) - 50.0) 10, 10, k 
CONTINUE ■ ^ ^ ^, 

GO, TO. 5 \ ' > ^ . 

W = P(I) 



We haver had to use the loop counter I- 
as an increasing counter and then use , 
I to define *J, which' is like a j 



/ 



decreasing .cOlinter, 



Ih. Label Statement 



X 10 



7 
20 



DO 10 r = 1. . ^ ' 

IF(ABSF(P(I)) - ABSF(TM)) L ip /aO 
IF(ABSP(pCi))^ . ABSF(T)) lO,^10, 5 
T = P(30 
QONTINUE 

iF(T) 7, 8, 7 . - . 

PRINT 20 

TOfeMAT(5HDN0NE) . 
STOP 



ERIC 



• 57 
■ 53 • 



15 « Label Stat ement 



10 

20 

3 



6 
7 

30 



DO 10 I = 1, N . 
IF (P(I) - ThV 10,%* 3, 
1.C0NTINUE J "^^^ 
PiRIIJT 20 . , . 
F0RMAT(5HDNdNE) 
STOP 

T = P(I) ^ ; 

IK = L.+ 1 • j 

DO 30 K = IK, N 

if(t - ,P(K)) 6> 30, 30' ' 

IF(P(K)*-'TM) 7, 30/30 

T = P(K) 

CONTINUE 



16, lAhel Statement 



k 
10 



SMALL = Q(L, 1) 
DO 10 J = 2, N 

IF( SMALL,- Q(L, J)) 10, iO, k 
SMALL =*Q(L, J) ^ 
CONTINUE 



17 . Label §tatjement 



10] 



DO 10 K" = 1, M 
I = M - K + 1 

iF(Q(i, IR) - T) ao, 3, 3 

CONTjte 

IROW'= 0' V.^ . 
GO TO 5 
IROW = I - 

BIG- = Q(I, IR) ^' * 



J 



ERLC 



\ 5 



Angwers to Exercises Fll-2 'Set C 



/ 



1(a) 



Kb) 



Label Statement or declaration 



50 
51 

10 
52 



50 
51 
52 



7 
10 

53 



DIMENSION X(50) 
READ 50, N 

f6rmat(I3) 

READ 51,, (X(I), I =-1, N) 
FORMAT(5F10.5) 
READ 51; A 
FNUM '»-, X(l) - A 
DO 10 J = 2, N 
mm = FNUM ^ (X(J) - A) ' 
CONTINUE 
PRINT 52> FNUM 
FORMAT(lHa,F15.5) 
STOP ' ' 
END 



Labey St^em'ent or declaration 




)NSION X(50) 

50, ^ 

tRMAT (13) 
-^51, (X(I). 1 = 1, N) 
^ORMAT (5FIO.5) 
READ 52, K, A ^ , 
FORMAT (13, FIO.5') 
FNUM = 1.0 
DO 10 J = 1, N 
IF(J - K) 7, 16, 7 
FNUM =: FNUM ^: (ic(.j) *^a)' 
CONTINUE \ ' " ^ 

PRINT -53, FNOM- 
F0RMAT(1Hd,F15.5) ^ . 
3T0P . . • 

END 



Label .Statement or gleclaration# 




50 
51 



^ 7 
10 

52 



DIMENSION X(50) 
READ 50/ N 
F0RMAT(I3), 

READ 51, (X(I), i = 1, N^ 

FORMAT(5F10.5) , ' 

READ 50, K 

DEN = 1.0 

DO 10 J = 1, N 

IF(J - K) 7, 10, 7 • 

DEN = DEN * .(X(J) X(k)) 

CONTINUE y 

J^R'ito 52, DEN 

FORMAT(1Hd,F15.5) 

STOP 

END 



5.9 



55 



3("b), ifiabel '< Declaration ^' 



j ; 




DIMENSION NP(1+) 


1 t 
1 




READ 50, (NP(I), I = 1, U) 

FORMAT( 1+115 ) 

READ 50, ICV , 

ISUM = 0 ' 

DO 10 L = 1,1000 

READ 50, NS, M 

K = M + NS - ((M + NS>/i+) * 1+ 
ISUM = ISUM + NP(K + 1) 


i ■ 50 

1 

/ ' 




io 




IF(XABSF(ISUM) - ICV) 10, 10, 8 




CONTINUE 


8 




PRINT 52, L, ICV, IsqM 
FORMAT( 1Hd,3I10) 
STOP 

END ■ , 
• 


. 52 





Comment:* When we gave the students a list of* FORTRAN function 
names in Table F2-2 we failed to- tell them about XABSF. 



v - 



r 




Answer io Exercise Fk-^ Bet A 

) 

Reyised flov chart 



1 




2 ' 


3 






"a" "b" "c" 




SHIP AIJD 




/ 














A, 9, 


C 


















h 












^B^ 





ERIC 



^^rogram ^ 

^^el Stat^ent or declaration 



A, B, C, D 




C 
C 



50 



10 
51 
3 
52 



93 



[INT 5C 

f ^Tl^^' iro, ^9x, iHC, ux, iro) 

A^JEPARATE PRINT ^^MENT IS USED !^PRINT A BLANK 

?RINT 51 1 ^W*" \ /' . . • 

FpRl# (iHo) ' ^ ' • 

A, B c 

F0RMAT(3F15.1) ^ • 
D = SQRTF(A*^#-B*^ + G^) 
PRINT 53, A,.B*/G, J). . . 
P0RMAT(1X, Fl5-i, ^2F10.1, FI5.I) 
GO TO 3 / 
END , . 



ft 



1 

^ Is 



Comment: . There Is another way to create the extra blank li^e,. It is 
with the use of the slash which is 'the end-of-rec.ord symbof^ed in ./ 
format code, k have decided to omit this in our presentati%of 
FORTRAN to^ students. Its use does not add enough to the^Under- v- 
standing of programming to warrant inclusiQrT--but you should at least 
be aware of its existence. If we want \o skip N lines after printing 
we simply add N slashes at the end of. the format. 

50- FOR14AT(iHl-,13X,lHA, ™ lHC,UX,lHD/^ckuses extra line to . 

be skipped 

In this way we eliminate the need for Staf ement 10 and its associated 
FORMAT declarations ' \ 

. 57 .0 t 



Ansvers to Exercises, Fk-^ Set A 
i. ' Label Statement /' 



5 
10 
20 



BIG = 0 

DO 20 I = 1, M 
IX) 10 J = 1^ N 

IF(ABSF(BIG) ABSF(P(I,J))) 5, 10, 10 
BIG = P(I,J) 

CONTINUE ^ " , 

CONTINUE ' . 
PRINT 50, BIG, 



Label Statement 







TLARpE = P(l, 1) 






JCROW = 1 






ICOL =1 

DO 20 I = 1, M 






DO 10 J = 1, N ' 
IF(TLARGE-P(I,J)) 5, 10, 10 






5 




TLARGE =P(lyJ) 
IROW = I'i 

ICOli = J ' 


10 




CONTINUE 


20' 




CONTINUE ^ 






PRINT 50, TLARG^ IROW, 'ICOL 



. 6 

7 

10 
20 



Label Statement 



ZLEAST = 0. ' . ' 

ZTALY =0. 

DO 20 I = 1, M, 2 

DO 10 J = 2, N, 2 t 

IF(P(I,J)) 6, 5, 6 

ZTALY za^ALY + 1. 

GO TO io ■ ! 

IF(ZLEAST - P(I,J)) 10, 10, 7 
ZLEAST = P(I,J)* 
CONTINUE / 
CONrmfE 

PRINT 50, ZLEAST, , ZTALY 



r 



Label Stater^ent 



% 




DO 20 I = 2*, M 






bp 10 J =• 1, N 

PCI,J) = P(I,J) + T * P(i; J) 






10 




CONTi^fUE^ 


' 20 




CONTINUE 



alternatively (and' better) 

DO 50 J - 1, N . 
TEMP V P(1,J)*T 
DO i^p I = 2, ^ 
P(l,j)-=: P(I,'J) '+ TEMP ' 

ifO Continue 

50 CONTINUE 



ERJC 



58 



62 



/ 



•5. Label Statement 



10 
20 



DO 20 J = 1, 

FMIN ='*P(l. J) • 

IROW = 1 

IX) 10 I = 2, M 

IFdMN ; P(I,J)) 10, 5, 5 

FMIN = P(I,J)' ^ 

IROW = I 

CONTINUE 

CONTINUE " 

PRINT. 50, PMIN, J, IROW 



6. L^bel Statement 



10 
^ ^ 20 



SUMl = 0.- ' 

IX) 20 I = 2, M 
L = I - 1 

IX) 10 J = i; L 

SUMl = SUMl + P(I,J) 

CONTINUE 

CONTINUE V 



0 



Cc-»nt,l: Students .re ll.M,.to b, o„al.„ on this proM» a. th. 
iteration box ^ ^ 



J <- 1 








J < I 




J <- J + 1 


w 



'is not the' same as 
a^d to say ^ 
is .illegal. 



DO 20 J = 1, I 
' . DO 20 j =^1, 

illegal 



Connnent 2: Ref4r ^o-tHe co~ in the Ite.chers , Connnent^ of the Main 
Itext on this problem. It is that |:ox;2'can be redrawn as ' 



i^i^ . 


I < M 


Qchanged from 2^ 


I + 1 


> 



^thout ch^ing the effect o|, t^e _ alglrithn. However, ve cannot make 
the corresponding change to the corresponding DO^tatement when using 
many of the FORTRAN II- processors. If wrote ' ' 



• DO 20 I = i,-,:m.-- 
then when entering the statement 



^ DO 10 J = 1, L . ' : . 

(wh^re L L - l) there will be nS initial test made to see if J which is 
i exc;eeds L which is 0. If the test were made all would be well. ' 

?• 1 Label Statement 







SUM^a = 0. 








MMiOC. = M - 1 




t 




DO 20 I = 1, MMAX 




1 

\ 




IPl = I + 1 




1 

i 




DO 10 J = IPl, M 






SUM 2 = SUM 2 + P(I,J) 




• ' 10 




CONTINUE 




. . 20 




CONTINUE • • 


r 


. 8. Label 


Statement* 





6 
10 

20 
10 

60 

■ 11 



lANY = 0 
LIMIT = M 2 
DO 20 L = 1, LIMIT 
J = M + 1 - L 
ZLAST =. P(1,J) 
JMr= J - 1 
DO id I = 2, JM1|# 
IF(P(I,/) -'2. * ZLAST). 6, -6, 7 
ZLAST = P(I, J) 
CONTINUE ' • . 

GO TO 20 • ^ 

PRINT 50, P(I> J), I, J- ' 
lANY =1 

CONTINUE . - ^ 
IF(IANY) II, 10, 11 
PRINT 60 
F0RMAT(5HaN0NE) 




» See the Teachers Commentary 'of the MaS^Text for alternate* f^ow 
solutions to Problem '8. « • ' ^ •• » > ^ 



chart 



...ERIC 



6o'-- 



'•,1. * 



Ansver^ to Exercise 'Sh^k Set 3 
. Label Statement or declaration 



C 
C 



ko 

50 
100 



NUMBER .OF NOIJCONGRUENT TRIANGLES WHciSE SIDES 
ARE OF INTEGER LENGTH LESS -THAN lOO ' 
IS = 0 ^' . ? ' 

iX) 50 I = 1, 100 
^ = 1 + l/2 
. K) ifO J ^ JJ, I ^ • 
IS = IS + 2*J - I 
CONa?INUfi 
CONTINUE ^ ? 

PRINT 100, IS 
FORMAT(llO) 
STOP 
END • 



(b) Label' ' Statement or declar ation 

ACCUMULATED PERIlfelERS FOR THE S TRIANGLES 



- C 

. C* 



30 
ho 
50 

lOG- 



COUNTED BY PRECEDING PROGRAM 
IP = 0 ' 
DO 50 I = 1, lOO 
JJ = 1 + l/2 
DO ifO J = JJ, I 

KK = I . J + 1 
DO 30 K = KK, J 

• IP = fp + I+J+K 
C0NTIJe5E 
CONTINUE . _ ^ 

' CONTINUE' 
PRINT 100, IP ' \ 
5X)RMAT(I10) 
STOP 
'END 



• -Answers to Exercises ?h^h Set C 



1, Label • Statement 



y 




READ 50, N 
FORMATf 15 ) 
FN = N 

LIMIT = SQRTF(FN) 










DO 10 K = 2. LIMIT 

IF(N - (n/K) * K) 10, \y 10 


3 








PRINT 51; K/ • 


51 




F0RMAT(lHa,I5) 




N =r N/K 

GO TO 3 


10 




.CONTINUE 

IF(N - 1) 7, 8, 7 


7 




PRINT 51^ N 


8 




STOP 
END 



* • '2.* Label* Statement or Declaration^* 



50 
5a 



^30 



1^0 
' 8 



10 
20 

52" 



DIMENSION A(^OO) 
READ 5©^N 
F0RMAT(I5) . 

READ 51, (A(I); 1 = 1; N) 

F0RMAT(1+F10.5) 

NMl = N - 1 

DO-20 J = 1, NMl . 

IF(A(J) - A(J + 1)) 20, 20, 5 

COPY = A( J) 

A(J) = A(J + 1) 

A('j + 1) =: pOPi i A 

tIMIT = J- 1 



IF(LIMJT) 30, 30, ^+0 
GO TO- 20 



DO 10 L = 1, LIMIT ' < 

K =' LIMIT +• 1- ^L 

IFCA(K) - A(K + 1))*20, 20, 8 

COPY = A(K) 

A(K) = A(K + 1) 

A(K + 1) = COPY 

CONTINUE 

CONTINUE 

PRINT 52, (A(I),'I = 1, N) 

F0RMAT(lHa,l+F10,5) 

STOP 

END 



prevents Inner DO from 
being executed when LIMIT- 



La"bel Statement 



50 
51 



10 
20 

52 



DIMENSION A(56o] 
BEAD 50, jfj 
F0RMAT(-I5 ) 
™ (A(K). K 
FORMAT .(J+|?ia^y 
NMi = N - 1 
IX) 20 I = -1, KMl 
MI = N - I 
DO 10 J = NMI 



1, N) 



1F(A(J) - A(J 
COPY = A(J) 
A(J) = A(J + 1) 
A(J + 1) = COPY 
tJONTINUE 
COOTINUE 
PRINT 52, (A(K), 
F0RMAT(J+F10.5) 
STOP 
END* 



+ 1\) 10, 10, 5 



K = 1, N) 



•4 



Label 
— 1^ 



50 
51 



^^2 

■ >6 
10 

:>i^ 
20 

52 



statement or declaration 



DIMENSION A(lOO), IB(lOO) 
READ 50; N 
FORMAT (15) 

READ 51, (A(I). I = 1, N) 
FORMAT (too;5) 
MAXINC = 1 " * ' 

DO ao J = 1, N' * 

B( J) ^= 1 ' * 

JMl = J - 1 ' 



iP(JMi) 36, 30, ho^\ 

GO ^ 13 / 



prevents Inner DO from being 
exeouted' when JMl = 0 



DO 10 K = 1 JMl 

iF(A(;t) . A(j)) 10, 10/5 

IF(IB(J) - llB\K) + 1)) 6, 10, 10 
IB(j)-= IB(K) +'3?^ . . . 
CONTINUE 

IF(MAXINC ^ IBCj)) Ih, 20>20 

MAXINC = B(J) 

CONTINUE « ' - . : . 

PRINT 52/MAXINC*'' 

FORMATS IHa,!^) . ' 

STOP ^ , \ ^ 

END ' * ' . 



SUBPROGRAMS 



, Use of the. computed GO TO statement 

k . 

■ FORTRAN Il-co^pilers for the smaUer coniputers somettoes do not allow US'' 
to define FUNCTION subprogrgps. All FORTRAN compilers, do, hovever, provi'de 
the so-called computed QO TO statement.* With this type of statement ve can 
simulate the.. action op^ functional reference as part of the main program. 
. We -give first a brief definition of .the computed GO TO and then use it' in 
simulating the FUNCTION subprogram given jLn Figure F5-2. 

Suppose the integer variable K is assigned a value, k, in the i-ange 
1 < k < m where m is some small number like 5". Now at certai'n points in 
■a program suppose ve would like to GO TO one of a set of m statements whose 
labels are s^, s^, s^. With each value k of. K we associate the ' 

statement label Sj^ ^ so that if k is Ihe". value of K, we GO TO label "' 
Sj^. The, FORTRAN statement that accomplishes this has the fom 



Example. 



00 TO (10; 19, 52, 6, 17), K 



Means: • Transfer, control to the statement labeled either 10, I9, 52, 6, or' 
17 when the Value^of K, is either 1, .2, 3, J^, .'or 5 respectively, 'if 
the value of K is not in the proper range, no gtiarantee is given as"\b » 
wBat will'be the "next task' executed. . . ' 

The flow cbiti^quivalent of this statement is- . ^ , 

i 



K< 0? 




K ^ 6? 



y ' ^ •/19 ■/ 



1 I 
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The computed GO "TO provides 'iJB with-e^convenie^^wltching, mechanism for 
tfeturning "6b a point of call, .simulating the automatic return 'that' is provided 
by the REfoRH -statement in a su"bprogram/* 

We now illustrate the switching mechanism for a program, much like the main 
•fl9w chart shotm. in Figure .^-'fi but 'which contains 'flfSjeddecH-n it the heart of ' 
tiie algorittmi used in the reTe^erlce flow cha^ in tj^.^same figure. ' There are 
two ^)Oints of call at, box 1 and at box k. Suy^>ose the detaij; of these -boxes 



Y ^ A + B X VX 



T 



, and 
^ I 
\ 



1 ^ 



P 4 X Vn 



The correspooding FORTRAN statements are 

Iy = A + E * SQRTF(X) ' 
'I r« ■ - ' 



T 



TEMP = N 

P =^4^. •^^^tfCtmp) ' 



* It 



Suppose the square root algorithm is" placecf* in the program beginning at 
the statement labeled 100. ^ Then we might' revise' statements. 1 andl h ^ as 



In place of "Statement l' 



51 



1 = 1- 
+ARG = X 
GO TO. IQO 
Y = A + "B * ROOT 



In pla(5^ of' Statement 4v / 
k 



61 



1 = 2 
ARG-= N 
CD TO 100 • 
P =» iv, ROOT 



The^ coding at statement 100, would end.in a computed GO TO as shown below: 



100 
2 

10 
5 



G = 1,0 ^ ^ ^ ^ ' 

H ^ 0.5 * (G +-ARG/G) • 
IF(ABSF(H - G) - ,0001) 5, 10, 10 
G = H - 

GO TO 2 ^ - " . 

GO 030 C51, 6l), I 



Retiirns control to Statement 51 
or 6ivdepending on the poiitt 
of call,' " ' 



. Answers ttf Exercises F3-1 



PUNCTION CUBRT(A) 
G = 1 

H = (2, * G + A/G**2)/3. 
IF(ABSF(H - G) - .OOOlO 5, k, k 
G = H 
GO;TO Z 
CUBRT = H 

RETURN * ■ " V - 

END' 



alternate 



FUNCTION CUBRT(a) 
G 1 

CUBRT=(2 ,*G+A/G*G)/3. 
IF( ABSF( CUBRT- G) - • 0001 ) 

' 5, 

G = CUBRT 

GO TO 2 . ^ 

RETURN 

END 



TF5 



FUNCTION FUNCT(X) 

FUNCT = (3. ^ X - 2.) ^ X + 1. 

RETURN 

END * I 



3. 



FUNCTION ABSOL(X) * 
IF (X) 2, 3, 3 
ABSOL = -X 
GfO^ TO 4 • 
ABSOL* - X 
RETURN 
END ^ 



Answers to Exercises F3-3 Set A 



1(a). 



'2. \ 



3U). 



\ 

^ C 
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FUNCTION ,FUNCT(X,Y) 

FTJNCT ^ ax^^3 +-Y)^^2 + 5- )/(ABSF(X) + 2.) 

RETURN • \ / ^ . ; 

END ^ - >y 



JJY^ FUNCT(R,S) + 6,,^ 



II 



FUNCTION RIGHT(A, B, -C) 
IF^C) 7, 7, 2 - 
IF (B) 7, 7, 3 ' 
IF-(A) k 

IF(C ^ C^-1a> a + b ^ B))V, 8, 5 
I?(A f A - (B ^ B + C ^ cn-g' 8 6 
IF(fe/ B (AV^\^ ^ M'7; 8 7 
RIGHT =0" , . ' V 

RETURN ^ 
RIGHT = 1 ^' 
RETURN 
END 

5TJNCTI0N XMAX (XJ Y, Z) ' 
XMAX = X * 
IF(XMA?C - Y*) 5, 3, 3 
IF/XMAX, - Z) 6, 4, 4 
RM33RN • ^ . ' 
Y 

GO 10 3' 

XMAT = * ' / . ^ 

RETORN 

END - * . 

READ 11/ A, B^C 
XLRGST = XMAX (A, 'B^ C) 
PRINT 11, XLRGST , 
STOP ^ 

FORMAa?-(3P15.8) 
END - y 

DMNITION'OF FUNCTION XMAX GOES HERE 
END 



f 
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70 



r 



!rF5 



If.*' Let a zero value for IQUAD represent 
coordinate axis. / ^ 



error exit for a point on a 



FUNCTION IQUAD(X, Y) 
IF(X),5, 2 
IF(Y)*.1|, 8,. 3 
IQUAD = 1 
RETURN 

IQUAD - 1+ • ' * 
RETURN 

IF(Y) 7, 8, 6 
IgUAD =2 ; , 

RETURN 
IQUAD = 3 
RETURN 
.IQUAD = 0 
RETURN - 
END 



6* 



r 
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a 

3 
5 

51 
6 

"61 

9 

^ 10 

%8 

11 

. '7 
12 



5 
10 
6 



FUNCTION INTSCT(X1,Y1,R1,X2,Y2,R2) 
IF(R1) 12, 12, 2 
IF(R2) 1^, 12. ^ 

DIST = SQRTfT(X2-X1)**2 + (Y2-Yl)**2) 

if(dist) 5, 10, 5 * * ' r 

IF(PIST - (R1+R2)) 51', 7, 51 

IF(DIST - ABSF(R1-R2)) 6, 7, 6 ; 

1F(DIST - ABSF(R1-R2)) 8, 6l, 6l 

IF(R1 + R2 - DISO?) 8, 9, ? 

INTSCT = 2 ^ 

RETURN ' • ' 

IF(R1-R2) 8, 11> 8 ^ 

INTSCT = 0 . 

RETURN 

INTSCT = iOO 

RETURN ' ; - ^ 

INTSCT = 1 . 
RETURN ' 

INTSCT = -1 ^ ' 

"RETUI^" ' ' ^ 

END 



FUNCTION ANYRT(Y,N) 
XN = N 
G = 1 

DO la I = 1, 10 

H = Y/G**(N - 1) ' 
-T^^ IF(ABSF(G - H)^ - .0001) 6, 5, 5 

G = ((XN - 1.) # 6 + hVxI^ 

CONTINUE 
ANYRT = G 
RETURN . / , 

END . * _ . 
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7 i 



10 



.FUNCTION TIRATE(N,R,TL)- 
TIRATE = TL 
RATE = -01 * R 
DO 10 I = 1, N. 

TIRATE = TIRATE * (l. 

CONTINUE 
RETURN 
END » 



RATE) 



•Alternate .solution: 



FUNCTION TIRATE (N,R,TL) 
TIRATE = ^ * (1. 4- R/100.)^* N 
RETURN 

END . ^ • * 

m 

Comment: ' * o 

- There i§ in FORTRAN II a so-called "statement ^nption" form' of 
^function definition vhich'vould be applicable to the alternate solution. ' 
(We do not feel it is important enough to inject this in the student text.) 
Using this statement- function form, 6ne can vrite the whole function 
definition as a single declaration: 

CaS^S"^? ^^^^^ DEFINITIOl^ 

jriRATE (N, R, TL) = fL *.(l. + r/iqo.) ^ N 



7(b), 



6 

' 7 
100 



9 
10 
101 
102 



ISWTCH =0 ' , 

= 100. / - * 

Y'= 100. ■ 
Dp 100 i = 1, 36 
X = X% 1.01 
Y = Y-+ 1.125 
IF (X-Y) 6, 6, 5 ' 
ISWTCH ='l 
MONTHS = I ^ 
IF ((1/12) * 12 - I) 100, 7, 100 
PRlto? 101, I, X, Y 
CONTINUE 
IF ((ISWTCH) - 1) 10, 9, m 
PRINT 102, MONTH ' 
STOP 

FORMAT ( 1HD,D15 , 2F1 0 . 2 ) ; 

FORMAT ( 10HaM0NTHSO=O, 12 ) 

END 



'5 
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Alternate solution 



6 
7 

'100 



9 
10 
101 
102 



ISWTCH = 0 
X = 100. 

Y = 100. •* • 

DO 100 I = 1^36 

X = TIRATE (I, 1., 100. ) 

Y = Y + I..I25 

If (X - Y^^ 6, 6,. 5 
ISWTCH = 1 ^ . 
' MOmjI = I 

IF((I/12) * 12 - I) 100, 7, 100 
PRIOT 101^ 1, X, Y 
CONTINUE ' . ^ . 

IF( ISWTCH - 1) 10,9^ 10 
PRINT 102, MONTH 
STOP 

FORMAT( IHD, 15 , 2F10 . 2) 
!, JX)RMAT( 10HdM0NTHSD= D, 12) 



7(c), Alternate solution 



3 
5 
6 
101 



X = 100. 
NX = 1 

X = X 1.01^ 
IF(X - 200.) 6, 5, 5 
PRINT 101, NX, X , 
STOP 

NX = NX + 1 ^; 
GO TO 3 

FORI4AT(5HdNXCJ=, I5,3HXa=,F10. 2) 
END 



Ansvers to Exercises F^-3 Set B 

1. * jLsibe:^ Statement or declaration 

C 
C 



. ERJC" 



5 



GREATEST COMMON DIVIS(^ AS AN INTEGER-VALUED 
FUNCTION CATJ.ro KGCD 
FUNCTION KGCD(KA, KB) . 
IF(KA - KB) 5^ 5, ^ 
KR KB / 

KB,= KAV - * 

KA = KR - • * . - 

jf(ka) 6, 6 

KR ^ KB ^ (KB/KA) KA 

KB'="KA' ^ - * / 

KA = KB • t 

GO TO 5' . • ^ ' ' % 

KGCD = KB^ , . . ^ 

RETURN g 

. ; 

70 



2. Label ^ Statement or declaration 



GREATESOl, COMMON FACaX)R ALGORITHM 
CALLED KGOP. IT CALLS ON KGCD 
FUNCTION KGCF(«A,KB,KC) * 
KX ='kGCl)(KA,KB> 

KGCD = KGCD(KX,K€)* ^ 

RETURN 

END. 



/ 



TF5 



3(a): 



Label Statement or declaration 



6 

ko 

50 
100 



NUMBER Oi* NON-SIMILAR TRIANGLES AMONG 
THOSE OF INTEGER LENGTH! LESS THAN 100 
IS = 0 N 
DO 50 I =-^*l,-100 
J J - 1 + 1/2 
DO kOJ = JJ, I 

KK = I - J + 1 

DO 30 K = KKy J ' 
IF(K33CF(I,J,K) ^- 1).30, 6, 30 
IS = IS + 1 ^ 
CONTINUE . 
CONTINUE , '^ 
CONTINUE 
PRINT 100, IS 
FORMAT(IIO) 
STOP ' 

ISND ^ . * 



3('b). Statement 6 to be replaced by 

6| Ins = NS + I + J + K 



Comment: In case students nm this .problem and also Problem h on 
the computer, you will probably want them to cut dovn the size of 
the pl*oblem,otherwise excessive computer time may be required. 
For exaa^le, triangles whose lengths are less than ,50 might be 
illustrative enough. Likewise in.Pro^em h, you could consider 
all numbers less than , 10^ instead of 10^. 



ERIC 
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LabeJ. Statement or declaration 



7 
8 

9 
10 

11 

12 
101 

13 



50 
60 



DIMENSION ICUBE (999) 
DO 60- 1 1^ 999 

ICUBE(I) = I 3 
DO 50 J = 1, I 

ITEST = ICUBE(I) + IOJBEC'J) 
IP(ITEST - 1000000000) 6, 60, -60 
K = I - 1 
L = J + 1 . 

IF(L . K) 8; 8, 50 ^ — 
.IF(ICUBE(L) + ICl©E(K) - ITEST) 10, 9, 11 
IF(KGCF(J, K, L) - 1) 13, 12, 13 
• L = L + 1 

GO TO 7 
K = K - 1 

GO TO 7' ^ • 

PRINT 101, I, J, ITEST, K, L 
FORMAT(l6,9HDajBEDD4-a, l6,9HDCUBEDD=a, 

110, 3Ho=a, l6, 9HaajBEIK:H-a I6, 6HaaJBED) 
L ^ Jj + 1 / 
K = K - 1 ' • ^ ' ^ 

•GO TO 7 ' 5 

CONTINUE '"'^ 
CONTINUE 

STOP 
END 



72 



TF5 



Answers to Exercises V3-k Set A 



r 



1. • 



2(a). 



SUBROUTINE ABSOli (X.ABSX) 

IF(X) 2, 3/3 

ABSX = -X 

RETURN 

ABSX = X 

^RETURN ' 

END 



SUBROUTINE CXADD (Al, Bl, A2, ^, A, B) 
A = <1 + A2 77,/ 

B = Bl + B2 

RET^ 

END 



SUBROUTINE CXSUBT (Al, Bl, A2, BS, A, B) 
A = Al - A2. ^ f f , J 

S = Bl - B2 „ ' 

RETURN 

END 



SUBROUTINE CXMULT>'(A1, Bl, A2, BS, A, B) 
A^= Al * A2 - B1-* . ' 

B' = Al * B2 + A2 ' ^ ' 

RETURN 

END . ^ . 

t 

SUBROUTINE CXDIV (A1, Bl, A2, BS, A, B) 

Jmm = A2*A2+B2*Be 

A = (Al * A2 + Bl * BS)/DEN0M 

B = * Bl - A1 * B2)/DENC»I 

RETURN * 

END ^ > ' 



«» 
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11 

5 



8 

9 
'12 



READ 11, Al, Blf A2, B2, OPER 

PRINT 11, Al, Bl, A2, B2, OPER 

FORMAT (5F10.I^) 

IF (OPER - 2.) 1^, 5, $ . 

CAL5, CXADD (Al, Bl, A2, B2, A, B) 

GO TO 9 

feALL CXSUBT (Al, Bl, A2, B2, A, B) 

GO TO 9 ^ ' ' ' 

IF(OEER 3.) 7, 7, 8 " 

CALL,CXMULT.(A1, Bl, A2, B2, A'. B) 

GO TO 9 ~ 

CALL CXDIV (Al, Bl, A2, B2, A, B) . 
PRINT 12, A, B 

FORMAT (1HOFLO,1^,1H'+ FLO.l^.lHl) 
GO TO 1 . ^ . • ^ 

END , " ' 

DEFINITION OF CXADD, .CXSUBT, CXMULT', CXDIV GO HERE 
END Q 

73 ^ 
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K < 0. 



3. A normal exit is indicated by ERROR = 0, If ERROR = 1, then K< 
SUBROUTINE SORT 2(K, A, B. ERROR) 



Ma). 



2 
3 

50 
5 



1 

1, LAST 
Ad+l)) 



50> 50, 5 



DIMENSION A(lOOO), B(<L000) 
IF(K) 6, 6, 2 
ERROR = Q 
LAST-= K 

DO 50 I = 
IF(A(I) - 
CONTINUE 
RETURN 
COPY = A(l) 
A(I) = A(I + 1,) 
A(I + 1) = COPY 
COPY = B(I) 
B(I) = B(I + 1) 
B(I + 1) = COPY 
GO TO 3 I 
ERROR = 1 
RETURN 
END 



SUBROUTINE COUNT (N^NOFAC) 
IF(N) 5, 5, 1 

NOFAC =0 . » 

XN = N * .V. 
IBOUND = SQRTF(p) - 1 
DO 3 K = 1, IBOUND 
IF(N - K * (N/K)) 3, 3 
h\ NOFAC = NOFAC + 2 ^ ^ 

3 CONTINUE . ■ ' ^ ; 

IF (N - K * K) 8, 7, ' , . 

7 NOFAC = NOFAC + 1 ^ ' " . 
RETURN 

5 NOFAC = -1 ^ 

RETURN \ 
END * * . , 

Note: Since the argument of SQPTP' must be floating point, we 
XN as the floating point representation of N, 



(b). 



3 

11 
2 



DO 2 N = 1,1000 > 
'CALL COUNT (N, IFAC) 
IF(IFAC^2) 2, 3, 2 
PRINT ITTn, 

FORMAT (no) , • 

CONTINTte 

STOP 

END 

DEF£NITI0N of ''SUBROUTINE COUNT GOES HERE 
END 
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5(a). 



(^). 



10 

11 



4. 



8 
10 

9 

11 
21 
22 



SUBROUTINE ALQUOT (NUMBER, N, IPART) ' 
DIMENSION* IPABT (50) 
rE.(NlJMBER) 8, 8, 2 
N = 1 

IPART (ly = 1 
IF (NUMBER - 3) 9, 9, 
XN = NUMBER, 
IBOUND = SQRTF(XN) 
DO 10 K = 2, IBOUND' 
IF (NUMBER - K * (NUMBER/K)) 10, 7, 10 
N = N + 2 
IPART(N - 1) = K 
IPART (N) = NUMBER/K. 
CONTINUE 

IF(P;ffiT(N) - PART(N - l))/9, 11, 9 
N = N - 1 
RETURN 
N = -1 
RETURN 
END 



DIMENSION IA(50) 
DO 10. I = 1, 500- • 
CALL ALQl$0(r (l, N, IA) 
IF(N) 9, 9, ^ 
ISUM = 0 • . 
DO 6 J. = 1, N 

ISUM = ISUM + IA(J) - ^ 

CONTINUE 

IF (I - ISUM) 10, 8, 10 . 
PRINT 21, I 

CONTINUE . ' 
GO TO ,11 
print' 22 
STOP 

FORMAT (110) 

FORMAT (llH IMPOSSIBIiE) 

END . • \ 

DEFINITION OF SUBROUTINE ALQUOT GOES HERE 
END . 




The first five perfect numbers are 6, 28, li-96, 8J.28, 33550336. 



■/ 
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- ■ .. . ■ \ -. . . . 

5(c). Let B be an array. B. is the sum of the aliquot parts for numbei* i/ 

' \ ' — ^ 

DIMENSION IA(50T, iB(500> ^ • 

DO 50 I = 1, 500 V ' " ^ ■ — 

CALL ALQUOT <I,N,IA) ^ * s 
IF(N)^1, 11, k ' f, 

ISUM = 0 

DO 6 J = 1, N . ^ : 

ISUM = ISUM + IA(J) 

CONTINUE ' a * 

IB(l) = ISU14 

IF(ISUM - l) 9, 50, 50 -' \' 

IF(IB(ISUM) - l) 50, 10, 50 ' * — 

PRINT 21,. ISUM, t 

^CONTINUE" * ^ ^ 

GO TO 12 

PRINT 22 . ' ^ -^—^ 

STOP . I ' '\ ^ 

FORMAT. (2110) . ^ ' 

FORMAT (IIH IMK)SSIBLE) 
END ^ 

DEFINITION OF SUBROUTINE ALQUO3? GOES HERE 
END ' J ' 



10 
50 



. 11 

12 



Alas, the only pair of friendly numbers, less than 5OO are^ 220 and 281+ 
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Answers to Exercises I^-k Set B 

1. d--- 



k 
60 



THE ^ LEAST FUNCTION IS CALLED TLEAST 
FUNCTION TLEAST (N,A) - 
DIMENSION A (100) 
S - A(l) 

DO 6o 1 = r. N 

IF(A(l) -^S)* 60, 60, h 

S = A(?) ^ 

CONTINUE ^ ^ 
TLEAST = S ' 
RETURN 
END 



60 



THE SUBLEAST FUNCTION IS GALLED MLEAST 

FUNCTION MLEAST (N, A) 

DIMENSION A(IOO) ' 

S = A(l) ' 

K = 1 

DO 60 = 1 N 

IF(A(I) - s) So, 60, h ^ 
S=.A(I) 
K = I 
CONTINUE 
ST = K 




END- 



3. C 



60 



THE PROCEDURE CALLED MARKS 
SUBROUTBIE MARKS '(N,A,S,K) 
DIMENSION A (100) 
S = A(l) 
K = 1 

DO 60 I = 1. N 

- IF(A(I) - S) 60, 60, k 
S = A(I) 

/ CONTINUE 
REirURN 
END 



Answers to Exercises F^-^ Set C 

1* Label Statement or declaration 



THE PROCEDURE CALLED DEGREE. 

FOR FORTRAN* S SAKE, SUBSCRIPTS ^ SHIFlED BY ONE. 
SUBROUTINE DEGREE- (N, lA) 

DIMENSION lA(lOO) ' " ' 

IF(IA(N+^) k, 2, k 
N = N « 1 " - ' ' 

IF(N)U, 1;- 1 
*RETURlf 
EHD \ 

80. 
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2. Lab^el ,Stabement oi^ declaration 



70 

60 

7 



PROCEDURE GALLED SMPLFY 10 

SIMPXiIFY THE COEIFFICIEOTS OF AN mm " 

DEGREE- POLYNOIvH/UL, N LESS TIIAN 100. 

(FOR FORTRAN'S SAKE, ALL SUBSCRIf^TS ARE 

SHIFTED OICS". ) ) 

SUBROUTINE SMPLFY (N,IA) ^ ( 

DIMENSION 4A(100) 

TES'f TO SEE IF DEGREE IS NEGATIVE. IF SO, DO NOTHING 
IF(N - (-1)) 7^7, 1 
ID = ABSF(IA(1)') 

NPLUSl = "N 1 * ^ • 

DO 70 I = 2, NPLUSl ^ 

ID = KGCD(ID, ABSF(IA(I))J 
' IF(lD - i) 70,^-7, 70 ^ ^ . ' 

COOTINUE 

DC5 60, I 1, NPLUSl ( . 

IA(I) = 'iaTi)/iD' ' 

CONTINUE f 
RETURN J . ^ 

END • , ' 



V*. 



Lat)el Statement or declapration 



C 
C 
C 
C 
C 
C 
C 



11 
2 



6 

60 
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PROCEDURE RDCMOS FOR REDUCING 
AN NTH DEGREE PpLblCKCAL, A(X), 
MODUH) AN MTH DEGREE POLYNOMIAL B(XJ 
WHERE M LESS THAN OR EQUAL TO N, . 
AND N LESS THAN 100. 

(FOR FORTRAN'S SAKE, ALL SUBSCRIPTS ARE 
SHIFTED BY ONE FROM THE VALUES IN THE 
FLOW CHART) 

SUBROUTINE.^ RDCMOD'' ( N,M, lA, IB) 
DIMENSION Ia(IOO), IB(IOO) 
IF(M) 10, 10, 2 ' 
IFCN-M) -10, 11, 11 ^ ^ 
NPl = N + 1 
MPl = M + 1 

IX = kgcd1;ia(M), IB(MPI)) 

IC = IB(MP1)/IX -> V 

ID = IA(NP1)/IX " 

DO 60 -I = 2, NPl 

FORM VALID SUBSCRIPT FOl^ N ■■ I 

NMINI = N - I 

IF(I - MPl) 5^ -5^ 6 ' 

FORM VALID SUBSCRIPT FCfe M ■■ I 



- I 



= IC 



MMINI = M 
IA{ NMINI). 
GO TO 60 
lA(NMINl) 
CONTINUE 
N = N - 1. . 
CALL DEGREE (N, IA) 
CALL SMPLFY (N, IA) 
GO TO 11 ' 
RETURN 
END 



IC ^ lA(NMINl) 
^ lA(ieGNl) 



ID^IB (MMINI ) 



ERLC 



Lal^ei Statement or declaration 



c 
"c 
c 
c 
c 
c 



11 

12 



111 

15 
16 

17 

ioi 

102 
"103 



PROGRAM FOR FINDING THE GREAT^ COMON 
DIVISOR OF ?W0 POLYNOMIALS * , 

THT'WO POLYNOMIALS lA AND 
IB HAVING DEGREES N AND 
M, RESPECTltTELY, AND EACH IS ' . 

LESS THAN 100. ^ - 

FOB- FORTRAN'S SAKE, INDEXES ARE SHIFTED BY ONE 
DIMENSION lA(llX)); IB(lOO) 
READ 101^ N, M 

NPl = N + 1 0 ' 

MPI = M + 1 ; ^ 

rlOl, (IA(I},, I = 1, NPl) ^ 
101, (IB(I), I =.1, MPl) 
CALL DEGREE(N, IA) I 
CAII/SMPLFY(N, IA/ 
CAI3?DEGREE(M> IB) 
CALL SMPLFY(M, IB) 
ISW = 0 ' ■ , 

GALL RDCM0D(N,M,IA,IB) - 

IT = N -\ ^ 

IFCjt) 15, ll+, 12 . ' ' ^ 

isw ^ 1 - "ISW ' 9 r 

IFdSW - 1) 7^'^, T " ' 

CALL TDCMOD (M,; N, IB, IA) 

IT = M " o ^ 

GO TO 11 ^ ^ 

^RINT' 102 

STOP ' ^ ^ ^ . /' * 

IF(ISW) 17, l6, 17 
PRINT i03, (B(I), I 
STOP ^ ^ 

PRINT 103, (A(I), I 
STOP 

F0RMAT(aDl6) 
J0RMAT(2H 1) 
F0RMAT(1H , 10l6) 

Eto 



1, MPl) 

^, NPl) 



Supplementajry Exercises in the Teacher's Commentary at the end of Section 

Here are the FORTRAN programs for the^fildw charts given in the solution. set 
to these exercises/, \^ ' 



1. 



.SUBROUTINE INTO NA, NB, NDEC) 
DIMENSION NA( 100) 
NDBC,=:. 0. :y 
DO^ 3 I = 1, N 

NI^EC = NDEC + NA(I) * NB**(I7 - l) 

'CONTINUE 

RETURN ' 

END 



PRir 



'19 
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' In*many implementations of FORTRAN there is no' parallel to the assign 

ment XOMP^-*-"0", Instead we must input the alphanumeric elements of 

vector NCCMP in the main program and include NCOMP in the par^eter ^ 

list fer function IDENT^ * ' - ^ 

SUBROUTINE IDENT (K, NA/ NCOMP, IVALUE) 
DIMENSION NA(lOO), NCCMP(l6) ^ 
DO 3 J^= 1, l6 , 
IF(NA(K) - NC(MP(<ff)) 3, H 3 
3 CONTINUE 
PR^NT 11, 

11 FORMAT (20HaINC0RRECTaCHARACTER) * ' 

STOP 

k IVALUE = J - 1 ^ - . _ 

RfcrURN. ' . - 

END - . , * 

, , A,'NC(M>, NDBG) J 

DIMENSION NA( 100)1, NeTMT ( l6) , , 

n6ec = 0^ ' . 

m'k 1 =,1, N . ' ' - 

^ALL nrar (i, na, nccm', nvalue) ' . 

NDEC = NDEC + NVALUE * 16**(N - I) <\ ^ ' 
RETURN 
END 



11 

12 



SUBROUTINE OUT (ND, NB, IR, M) 
DIMENSION IR(IOO) 

NQ = ND/N9 

IR(M) = ND - NQ * NB 

IF(NQ) 5, "5^ fe 

M ='M + 1' _ ' 

ND ^m'' 

GO TO 2 ^ 
RETURN 



DitalNSION NA(lOO), Nr!(100) 
READ 11, NBl, NB2, N, (NA(K), K = 1, N) 
CALL" INTO (N, NAy NB1,.NBAS10) 
CALL OUT (NBASIO j NB2, .NRl M) 
REVERSE ORDER Of1nRvENTRI|1S IN NEXT SIX STEPS 
MTODLE = (M + i)/2 
DO 10 

k = M + 1 ^"l 
NEED = NR(I) 
NR{I) = NR(K) 

NR(K) = NEEt) • , t ' 

CONTINUE 

PRINT 12^ (NR(I), I = 1, M) 
W TO 1 » ^ ^ 

FORMAT (316, (7151)) • 
FORMAT (iHnlOOIl), • , 
END . ' ' 

DEFINITION I OF INTO AND OUT GO .HERE 

END' : 



I 1^ MIDDLE 



7 ^ 



TF5 



6(a), As in the hexadeciraal-to-aecimal-conyersion exercise (Problem 2), the 
alphanumeric coraparison' vector, in this case ROMAN, intist b.e input into 
^ the main program and then transferred to the subroutine subprogram as 
one^of the components of the parameter list. Input the elements of / 
ROMAN in this '^order : I, V, X, L, C, D, M. • ^ * 

SUBROUTINE RNUM (N^A^.RO^, NUM) . V 
DIMENSION A(20), ROIAlJ(7), NVALUE(7) 

NVALUE(l) = 1 - ^ ^ . ' 

NVALUE(2) = 5 - . 

NVALUE(3) =10 . . * 

NVALUE(lf) = 50 

NmUE(5) =100 . . _ 

NVALUE(-6) = 500 ^ 

nvalue(7) = loot) 

NUM = 0 
LAST =8 

DO 50 K = 1, N . ' ^ - ' 

DO 6 I = 1, 7 ; ' - 

IF(A(K) - R0MAN(l))6, 7, 6 
CONTINUE ' o ■ 

PRINT 21 ' 
STOP , ' 

IF(last - I) 8, 9v. 9 , 
NUM = NUM - 2 * NVALUE(LAST) ^+ NVALUE(i) 
GO TO 10 ■ ^ , ^ 

NUM V NUM + NVALUE(I) * ^ ' 

LAST = I ' ' 

CONTINUE . ° ' - \ 

RETURN , ' - 

FORMAT • (20HDlNC0RRECItjCHARACTER) ' , 
END • * 



(b). 



•7 

8 

9 

10 
50 

21 



. 1 



..21 
22 

25 



DIMENSION A(20), B(ao), RmAN(7) 
READ 25, (R0MAN(I>, I = 1, 7) ' ' ' 
READ il„N, (A(K). K = 1, N) 
MpRINT 21, N, CA(K), K = 1, N) ' 
READ 21, M, (B(K)- K = 1, M) - 
PRINT 21, (B(K), K = X, M) 
CALL RNUM (N, A, NUMl) 
CALL RNIM (M,.B, NUM2) 
SUM*= NUMl + NUM2^ 
PRINT 22, SUM '-| , 

GO TO 1 f 
FORMAT (aO, 20A1) 
FORMAT (7HDSUMD=DFlG:b) 
FORMAT (7A1) 

END . ^ I . ' ' 

DEFINITION OF SUBROUTINE RJJUM QOES^HERE 

END - ' 



ERIC 



Si 
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Supplementary remarks on assignment of values to non-local variables 

^ The following should'^lielp you to visualize how ^the information on non- 

local variables is transmitted to a FORTRAN subprogram. Suppose a calling 
argument is a single v^iable, like T, and suppose it Is matched with a 
dunHny variable, say X in the subprogram. ^ A1j the time the subprogram-is 
executed, the subprogi^am "takes note of" the address where the value of T* 
ma^^be found (i.e., location of the window box labeled T). If the calling 
argument, is an expression like 2. * T, ' the subprogt'am wi^l instead take note 
of the address where the value of 3. * T has been temporarily stored. The 
programmer is normally not aware that the compiler, in generating the code 
for a call ^n a subprogram, generates instructions to evaluate argument'Vx- 
press ions and to save their values in "temporary" storage locations. Althoiagh 
the programmer is not told/whete these "temporaries" are located, their 
addresses are nevef the less made available to the subprogram each^ time the 
subprogram is called^. ^' . ^ ^ ' 

, We can now begin to perceive how the foregoing explanation might relate 

to our sealed chamber models with slips of papery boxes beipg dropped into 
the funnel. If a subprogram is told where the value of an argument is located, 
it can alter the value of the argument. Hence this corresponds in our model 
to dropping into the funnel the window box that contains the value of the , 
argument. If the argument is 2. * T, however, the subprogram is not told ' 
where T is located, only where the Xalue of ^ computation 2. * T is located. 
Hence the subprogram cannot alter the value of T.- In a way this siaggests 
the protection which is inherent in the slip of paper concept. Strict^/ 
speaking we could "protect" any variable used. as a calling argument to a* 
FORTRAN function by first assigning the , value of that variable to an auxiliary 
variable aijd then usijig thfet auxiliary variable as the argument. In this way 
no matter what the function did with the window box it received, it could not, 
."harm" the value of the original variable. 

Example ; Suppose the FUNCTION is called FUN .and it has arguments A, B, C. 
Then on^ way to simulate the slip of paper concept of protection is by writing 




instead of 



/ 3 



' Z ='fUN(T1,.T2, T3) * 
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A 'more sensible way to accomplish the same objective (and actually a 
safer and more efficient vay) is to let the subprogram itself reassign in- 
coming values 'of the dummy variables to new auxiliary variables that are 
strictly local'to the ^ subprogram. After, .this reassignment the auxiliary 
variables are used in place of the dummy arguments. Thus for the above . 

example, thfe call to FUN could be safely written, a3: 
# 

Z = FUN (A, B, C) * * • 

pi-ovided the definition of FUN is coded something like 

FUNCTION FUN (DUMMYl, DUMMY2, DUl^lMYB) 
ARGl = DUMMYl 
ARG2 = DUI^Y2 
AKG3 = DUMMY3 



0 



^gtion name arguments in FQKTMN II processors 

Some of the processors for the larger computers like the IBM 7090 do 
in fact provide an "improvised" way to define f\mctions having function name 
arguments (but no^ label arguments). To use a functior]^ name' as an argument in 
such procej^sors, it is. necessary to prpvi'de a special card in the calling 
program .which" lists the function names that 'are -to be used as arguments. This 
card has the letter F punched in Column 1. " ' 

There are other details' regarding the dropping of the. teiminal letter F 
-on librarjJ subroutines that are used as arguments. SucL details ile messy. 
Even if'ycfu have this capability available to you in the laborator^ one opuld 
question the, wisdom of .teaching it /o ^students at this level. In more advanced 
languages function name arguments are -handled in a more natural vay, ^'or fee 
details you should, course, consult the appropriate reference manual<^ 





5^ 



SUBROUTINE -FOR TWO SIMULTANEOUS EQUATIONS 
SUBROUTINE R00TS2(A1,B1,C1,A2jB2,C2,X1,X2,N) ■ 
DENC^ = Al * B2^- A2 * Bl ' 

IF(DENOM) 3, 5, 3 * 
XI = (iBi * C2 - B2 * Ca)/DENOM . « • 

X2 = (Al * C2 - A2 * C1)/DEN0M J 
N = 0 • ' ^ 

RETURN 
N = 1 
RETURN 
END 



3 



10 



11 
12 



13 
Ik 



A SUBROUTINE SUBPROGRAM FOR REAL 
ROOTS OF QUADRATIC EQUATIONS . 
SUBROUTINE ROOTg (A, B, C, N, XI, X2) 
IF(a5 2, 11, 2 
IF(B) 3, 8, 3' 
DISC = B*B-iv. *A*C 
-IF (DISC) 10, 7, 5 
N = 2 ■ i \ 

XI = (-B+SQRTF(DISC))/(2, * A)/ \ 
X^ = (-B-SQRTFCDISC))/(2. * A) 
GO TO^lk^. 
N = 1 

XI = -b/a, 

GO TO llf 

N = 3 > ' 

GO TO ll^ * ' 
IF(C/a).9) 9, 10 ' 
fJ = 2 

XI = sqrtf(-c/a) 

,X2 = -XI . jf 
GO TO Ik . j : \ 

IF(B) 12, 13, 12 I 
N = 1 * ' , . 
XI = 7C/B 

GO TO ll^ ' ^ ' 

N = 0 . 

RETURN,. 

END * ' * ' . 



A PROGRAM FOR ROOTS OF QUADRATIC EQUATIONS 
READ 20, A, 0 " 
raiNT 20, A, B, 0 j . 

CXlL roots (A,B,C,N,X1,X2) 
- 1) 5, 6, 11 

' 2) 99, 7, 8- - * . 

PRINT 21 
GO TO 1 
PRINT 22, 'XI 

GO TO 1 ' ' 

PRINT 23, XI, X2 

GO TO 1 

PRINT 2^ 

GO TO 1 

^P 

FORMAT(3Fl6.8) 

F0RMATf2i^HI3N0DlNTERESTINGaS0LUTI0N) 
FORMAT(l7Htl0NEDSOLUnONbx=a,Fl6.8) . 

FORMAT(l9ICTWa3SOLUTIONSQXl=a,Flg,8,5HdX2==a,Fl6,8) 
FORMAT( 22B3SOLUnONSDAREaC(M'LEX) 
END ' . . 



k 
2 



SUBROUTINE SUBF(X, Y,T,K) 

IF(X - 2)2, 4, 2 i 

RETURN 

^^=^X 3 + Y)| 2 + 5. )!/(X « 2. )' 

'^END . I 



12 

50 
13 



CALL SUBF(R, S, Z, iCEEST) 
IF(KTEST) IS, 13, 13 
PRINT 50 
GO TO 60 

FORMAT( 21B3ZDCANN0T[3BaCQMPUTED) 
z = Zj 6. y $ ' I 
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Answers to Exercises 

4 



1. 



SUBROUTINE C0NTCH(N, IS, IC, ICOUOT) 
DIMENSION IS(lOO) 
ICOUNT = 0 _ ' 
M = 1- 

CALL CHEKCH (N. IS, M, IC, LOC) 

IF(LOC) k, 6, k 

ICOUNT =-ICOyNT + 1 • 

M = LOC + 1 

GO TO 2 

RETURN 

END 



The symbols ) and ( must be input for the main program and be 
included in the subroutine parameter list, say call them NPARl and 
NPAR2, respectively* Also let ERROR equal^ 0 for a correctly 
written expression, 1 for a negative counter value and 2 for ^ 
nonzero Qounter value at the end of the scan. 



ERLC 



6 

1 
8 
-50 

10 

il 



SUBROUTINE PAREN^ TN,* IS, NPARl, ,NPAR2, ERROR) 

DIMENSION IS(lOO) ^ . 

ICOUNT =0 

DO 50 I = 1, N 

IF(IS(I) - NPARl) 7, 7 

ICOUNT = ICOUNT -^1 ^ 

IF( ICOUNT) 6, 50, 50 

ERROR =1 ' , 

RETURN 

IF(IS(I).- NPAR2)^50, 8, 50 
ICOUNT = ICOUNT +1 
CONTINUE. • 

IFCJCOUNT) 11; 10, U 
ERROR = 0 ' 
RETURN 

ERROR = 2 . . - 
R^3!UBN 

!end 

SUBROUTINE C0NTST(N, IS, K, IC, ICOUNT) 
DIMENSION IS(lod), IC(lOO) 
ICOUNT = 0 
M = 1 

CALL CHEKBT (N*, IS. 'M, K, IC, LOC) 
IF(lX)fc) 6, . 
ICOUNT = ICOUNT + 1 
M ^ LOC + L I 
TO 2 




89 

86' 



r 
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SUBROUTINE AVER(M.N,IA, V) 

DIMENSION lA(lOOO) 

NUM = N - M + 1 

ISUM = 0 ' 

IF(NUM) 3/ 6, 3 

DO 1+ I = M, N 'V * 

ISUM =,ISUM + IA'Ci) . 

CONTINUE 

V = isum/num 
retOrn 

V = -50 

RETORN 
END . ' 



11 

12 



DIMENSION lA(lOOO) ^ 
READ.11,K,{IA(I),T = 1,K} 
PEAD 12,M,N 

CALL AVER (M^N^IA, AVERAGE) 
PRINT 12,M,N^AVERAGE 
GO TO 1 

F0RMAT(I5(114.I5)) 

F0RMAT(2IlQ,F10-2) 
END 

DEFINITION OF AV^R GOES HERE 
END ^ • 




Chapter TF7 " 
SaME ^^ATHEMA^PICAL APPLICATIONS 



Th^ material here will , parallel closely that in the student's Chapter ^7, 
Since no new FORTRAN .concepts have been introduced, discussions are limited to 
specific points regarding the exercises* ♦ 



1. 



C. 
-C' 

c' 



33 
1 

2 

3 

if 



1 



PRINT 1 ' ' , ' ' ^ / ' 

F0RliAT(5X/ IHA, 9X, IHB, 7X, 1+HEPSI, 15X, IfHROOT) 
CALL ZER0K(0., '2., .1, l) 
CALL ZEROK(.l, 1., .15,2) 
CALL ZEROK(0., 2., 3) 
CALL ZER0K(6., 2., .1, k) 
CALL ZER0K(0.1+, h.) .1, 5) i 
CALL ZEROK(0., 2., l.'^-h, 3) 
STOP 
END 

SUBROUTINE ZEROK OOES HERE 
THE RJNCnbN SyBPROGRAM CALLED RJNCT 
IS GIVEN^ JffiRE- \ 
JUNCTION RJNCT (X,K) 
IF (K-2) 1, 2, 33 
IF (K-l|) 3, 5 ^ 
FQNCT =^'(X * X^- 1.) * X - 1, . * 
RETURN ' * ' 

RJNCT = X + LOGF(X) ' 
RETORN 

luNCT = 5- - X - s. * simix) ' 

REinjRN 

FUNCT = (X * X - 3') * X - 2. 
RETURN 

FUNCT = ((X + 2.) X - 13*) * X + 10. 
RETURN 

BND X , . ■ 



Calculates results 



(1.) (b) ^1.3K 



Qy EPSI = O.L 



(20 ^0.60625, :2PSi-= 0.15 r ' j ^ ^ 
- (30 0.87500/ EPSI O.k; 0.9if565y EPSI ^ Ip" 
;/(if.) (Not availafc^le this edition) , , • . . ^ 



ib»)' Method is inatrbli cable. 
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DO 1 I«l»4a ^ 
X« I' 

Yl« Fll(X) 
Y2« F12CX) 
Y3« FKX) 

Y4« F2(X) \ 
Y5« F3(X) V 
Y6« ,F13(X1 - ^ ■ 
Y7« F4(X) 
?8« F5(X) 

PRINT 2»Xf^Yl»Y2»Y3#Y4fY5fY6»Y7fY8 
2 FORMAK 9F10«2) 
STOP . • 
END' 

f^UNCTION FllCX) 
Fll« (X**2 -2«)»X -5^ 
RETURN 

END , 

FUNCTION F12<X),. 

F12 • (IX+3«)*X-2« )»X-4« i 

RETURN 

^ND 

FUNCTION FKX) ^ ^ 

Fl« (C3«»X-2«)«X*«2*f'7«)»X-4« 

RETURN' 
ENCr^ 

•FUNCTION F2CX» ^ 
F2« (X»»2-1«)*X-1« 
RETURN 

END > 
FUNCTION T3(.X) . 

F3« (X-3«)*X-4«#(SINF(X))»*2 , 
/ RETURN 
END 

FUNCTION F13rX) 
F13« X,-S,INFCX)ycOSF(X) 
RETURN 
END 

FUNCTION F4<X) ^ 
F^«. X^LOGFtX) \ 
RETURN ) 
END • 1 

F'UNCTION F&(X) ^ ' ' 
F5« 5«-X-5#*SINF<X) 
RETURN*' 
. END _ 
END 



F. 



-3. *^(af 15.03 (or 15 to' the nearest .foot) 
(b) 15.6$ fe§t ' ' 



TP? 



C 
C 



2 
3 

101 

i 
102 



MAIN 

CALL ZEROKLd, LABEL, 0., 3.1^159/2., 1.E.1|,A) 
IF (LABEL) if, 2, if ^ ; J ' m 

CALL ZEROKL (2, LABEL, 0., 34^159/2., 1;e5^,B) 
IF (label) 1+,- 3, 1+ . ^ 
PRINT 101, A, B ' ^ ' . 
FORMAT (l7HaROOiaOFaF(X)DISD,B20.8, ' 
17HDR00TCX)PM(X)DISD, G20 . 8) - 

STOP » . ; 

PRINT 102 : . ' ,^ ^ 

FORMAT (27HDMETH0I)aisalNAPPLICABtEaP0R, 

15HDF(X),D0R,DG(X) ) 
STOP » ^ . . . ' 

END • ' ' ■ . 

' SUBROUTINE ZTOOKL G0B5 HERE 
FUNCTION DEFINITION FOLLOWS . . 
FUNCTION PUNCT(X,K) 
IF (K-X) 1, 1, 2 ' 
PUNQT = SINF(X) - .66666667 * X ^ 
RETURN' 

FUNCT = SINF(X)/COSF(X) - .1 *'X 

RETURN 

ENB - 




8 
101 

7 
102 



33. 

53 
1 

2 
3 

iT 

6 



- MAIN ^ ^ - . z 

CALL ZEROKLd/ LABEL, p., .5^ l.E^lf, Fl) 
CALL ZER0KL'(2, LABEL, .5, 1., l.E-1^, F2) 
CALL ZER0KL(3, LABEL^ .^07, 1, l.E-if, RG2) • 
CALL ZEROiCL(i+, LABEL, RG2, 1., I.E-lf, RGjj ^ 
CALL ZER0KL(5, LABEL, i^G3, 1., ^.E-lf, RGi|'.)-. 
CALL ZER0KL(6-, LABEL, -RGif, 1., l.'E-lf, RG5y-\ 

IF (LABEL) -7T-8r"7-^^ ■ - : 

PRINT 101 ^ ' ' 

FORMAT^ (5Ha00PS) , , 

STOP . ' ' ' , 

PRINT 102, Fl, F2, BG2, RG3, RGi^; RG5 



FORMAT f2E20:8/liE20.8) 
STOP f 
END ' ' 

^ffi FUNCTION FUNCT(X,K) 
?1«JCTI0N PUNCT(X,K) ' 
IF (K^2) 1, a, 33 . ' ' 
^IF(K-ii) 3,^, 53 
^IF (K-6) 5, 6- 6 
FlJNCT ^ X * SQRTFd - X''* X) - ^25 
RETURN 

HmCT ±: -SQRTFfl - X * X) - X' 
RETURN ' . * ; 

FUNCT = SQRTF(1 ^ X * X)^ - X 2 ^ 
RETimN ' ' ' ' ^ I 

FUNCT = SQRTFd fX^X) - X^«^3' 
R&IURN ' ; 

pyNCT,= SQRTFd- - X *-X) - -X ^ if 1 




SQRTFd X X) ^- X p \ 

V3 



9X 



Si 



* * • 

Axis^ts to Exercises ¥J^2 • ^ * 

.(a)"*"lf EPSI is too small for the accuracy of the machine being used, 
^ it may not terminate ' ' ' 

(b) If N is greater tlian 100, - -the storage set aside for T will 
be exceeded- ^ 

(c) Replace statement 8 with'' 

8 IF (N --100) 3, 7, 7 / ' ' 

► The remainder of the program could be * 

3 N = N + 1 

GO TO 2 , . . • 

* ' 7 PRINT 103 

103 FORMAT (30H ERROR TOLERANCE EXCEEDED) 
, 9 AREA = T(N) / 

' , ' PRINT 106; AREA . ^ • . 

106. FORMAT (8HnAREAn=D, E20.6) - ^ 

. STOP ' ^ ' 

• ^ END 

2- We may eliminate the calculation of 2^ for each n. Also, since only 
the purrent sum and the previous sum are needed, we may denote them bV 

AR^ a nd OL AREA^ respectively,, thus eliminating the use of subscrip-J^ed 

variables.. ^ 



r READ 

If FORMAT <E 20#8) 

• - - OLAREA = 0%5 **(X»0*+ OtST' 
H d ll - ' 

H « 1 • 
N ■ 1 

2 H ■ 2«H 
L ■ M-1 " s 
H - H/2# ' i' ■> 

' S ■ 0 - : 

^ FK - >: ^ 

3 S ■ S-M^O/dtO+FK^H) 

^ - AREA ^0.5^0LAREA + H*S 

; - IF ( ABSF{AREA-OriAREA)-EPSl)l+,5/5 

. 5 1F<N-100). 9#7#7 r . ' 

9 N » H41 ' 

• TZERO ".JiREA^ ^ • • ^ 

• GO TO 2 ' . 

* 7 PRINT 8 • . ^ 

8 FORMAT <30H IERR0R> TOLERANCE ^01 SATISFIEOi 

4 RRINT 6#EPSI#TX/^ ^ 
6 -FORMAT (2E20#8I( ' 

'/ STOP . ' ' 

. ■ ENO^ • 



For ^I = 6,$1, TL = 0,691+12. 
- For ^1=^0:001, Tl = 0.69339. 

erIc - . .» 



3. 

« 



Let LIMIT the maximum number of iterations to .be, per formed. 

• y * ' 

5 READ 1, LlfcD-. ' 
- 1 FORMAT {Tiaf . • ' ' ' . * 

PpiOT 10^ ^ 
10 FORMAT (12X,2HT1 17X,1+HT)TFF 16X,5HLIMIT) . \ 
TZERO = 0.5 * (l.O + 0.5) . ' 
• - N = 1 

2 M = 2 * M 
L = M -.1 
? = H/2. 
S = 0 , 

D03K=1, L, 2^ « 
.FK = K 

3 S = S +-l./(l.+pk* H) 
Tl = 0. 5 TZERO + H S 

, DIFF = Tl - TZERO 
IF (N - LIMIT) 9, h 
9 N = N + 1 

TZERO = Tl ' ' 
GO-TO 2 

if. PRINT 6, Tl, DIFF. LIMIT ' a 

6 P0RMAT(2E2O.8,I15) 

ST^- . 'a' - • 

END ' ) 

For LIMIT = '15, Tl = Ot693ll+. ^ 




h. Be/ore the ENP statement we may insert' GO TO -8, 



5. * 



PRIHT 10 
10 FORMAT (10X,1|HAREA 16X,1HN) 
^ h READ 1, n' 

1 FORMAT (110) j 

S '= (1. + l./2.^)/2. 
KLAST =-N - 1 
*XN = N 
DO 2 K = 1, ^JCLAST 
XK = K » 

2 S = S + l./(l. + Xig^XN) 
AREA = S/XN 

PRINT 3, AREA, N 

3 F0RMAT(E20.8,I15) ^ 
^ GO TO 1+ 

END 



For N. 5, I AREA 

For N ='^250 AREA 

For N = 75J AREA 

^*or N =^125, AREA 

For N ^ 206, AREA 




0.69562 
0.69325 
O.693IF 
0.6931 
0,69315 



95 
93 



Area under curve : the g^eneral cas-e 
Exercises F?*-3v ^ 



Answer^ 



1. 



(a) 



2«AREA2(0»fPn5000) 
PRINT 4t2 
4 FORMAT (IXEi 
STOP 
END 

PyNCTION AREA2(Af\B»N) ^ 
FN=N • ^ 
H=(R-A)/FN\ V 
S:^ (FUNCTIB)-4'FUNCT(A))/2# 

IF it) ?»2>3' 
3 DO 1 K=i»L 

FK= K . , 

1 S ^ S+FUNCT(A+FK*tf) 

2 are;^2= h*s 

RETURN ' t 

END ; 

, FUNCTION FUNCttX) 

FUNCT* SIMP(X) ^ 
c . RETURN I 
' END 

end: 



■ i 



The area under the sine curve Is 2.000"^ while for, a, semicircle of 
diameter Jt the area is 3.876. . ^ 



(b) I Of course this is not th^ usual method for pfinting.a logarithm 

lablef The prpblem is given to connect again the area method to 

'I •. * 

he introductory discussion of i,hx. The correct values fo: 



ire these: 






inx 


1 


0.00000 




1.T9176 


.3-1 ' 


2.^790 


l6 


. 2.fe9 ^ 


... 21 -V- 


3.01^^52' 


• 26 


3.25810 



X 

36 

11 



r t 



3.5535?\ 
3.71357 ) 
3.'eSB6U 
3.93183 ♦ 



2. The variable M gives the current number of subAivxslons . These 
, revisions would accomplish the desired termination. 

^ (1) * N, the-maximum number of subdivisions, may te included as^an 

argument of the .function:' AREA (A, B, EPSI, N) . • ^ 

f « 

(2) The>emainde;r^bf the program after Statement 5 could b^ 
*■ • ' 

. ' . - IF(M N) 6, 7, 7 

6 GO TO 1 V . * 

/ . 7 PRINT 8, Tl ■ - I 

8 FORMAT (7HnAREAD= E20. 8,30HDcxiACcWci^aCRITERIONaEXCEEDED) 
' 2 AREA = Tl ^ ' " 

RfiFORNt . ^ 
EM) . ' ' 



(a) 



00 4 LWi,3 . , 

Z«AREA2M#f3»»N) 

4 PRINT 3fNfZ * 

.3 F0RMAT<ikH0N«Dl3», 9HdddAREA«dE20«8) 
Z«AREA(l«f3««0«001) 
PRINT 5fZ 

5 FO1?MAT(6HDAREA«E20«8). / ' 

. - STO?^ ^ , 

. / . -^NO . . ^ ... • 

C FiJNCTIOH AREA2 GOES HERE ^ 
C FUNCTION^ AREA GOES HERE 
FUNCTION FUNCTUS* 
FUNCT«»43429/>«r 
RETURN. ; ' 

;eno - \ \ 
END 



Galcuiated results: 

Subdivisions EPSI ' Area 

i . " ' 0.57905 

? ^ . 0^50667 

10 ' 0. 1^7721^ 
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3. . (b) 



. C 
C 



DO 6 L2 = 1, 3 . 

N = 2^*(L2 - 

Z = AREA2(-2., 2., N) 
6 PRINT 3, N/Z * V 
3 F0RMAT(1+H N = I3, 9H " AREA = E20\8) 

Z/= AREA(-2,, 2., 0.001) 

print' 5, z . " ' 

5 "fORMAT(6H ^AREA = E20.8) ^, . 

STOP . ; . ; ' 

FyNC??ION AREA2. GOES HERE 
FUNCTION AREA GOES HERE ' ^ 

FCNCTION FUNCT(X) . . ' 
KJNCT = (3.^X + 2.) ^ X + 1. 
' RETURN 



END 
END 



Calculated results; 
. Subdivisioas 
1 

• ' 2 
1+ 



EPSI' 



10 



-3^ 



Area 

52.00 
28.OQQ 
22.000 
20,000 



/ 



/ 




(c) 

• \ ■ 



1 



C 

c 



Z«AREA2(1«»4«»N) 
7 PRIHT 3fHfI^ ^ 

FORMATt^H N« 13» 9H AREA« E20.8) 
2«AREAC1.»4»»0.00.1) , 
PRINT 5fZ • 
5 FORMAT (6HDAREA«g20«8) 

STOP . ' 

END 

FUNCTION AREA GOES HERE ^ 
FUNCTION AREA^GOES HERE 
FUNCTION FUNCT(X> 
FUNCT«(X-1'.)»X»X^" 
RETURN 



END 
END 



Calculated results: 
Subdivisit)ns 

. 1- 

2 

h 



EPSI 



10 



r3 



Are^ *\ 
•'72.00 

'50.063 
If u. 578 
1^,750 



ERIC 
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TF7-lf Simultaneous linear - equations : dev^mnpT ^^ systematic mejho^ ^ 
of solution 



/"^ Answers to Exercises g?-^ 

1. 



3 1 READ l/AU,AlX.Bl, A21, A22,^B2 
/ fRim 1. A11/A12, Bl, A21, A22, B2 

1 FORMAT"flXP9.0,5FlO.O) ^ - 
AAli ^ An/All 

•^AAia A12/A11^ 

"bbiXbi/aii 

A21 - A21 AAli 
^2 •= A22 - 'A21'-^ AA12 • 
Bfi2 ^ B2 - A21 * BEL 
X2 = BB2/AA22 
Xi = BBl - AAX2 * X2 
PRINT 2, Xl^ X2 

2 - FORMAT(7HCEDX7a=E20.8, 5Xl+HX2a= <E20,8) 

GO TO 3 




END ^ 

1.6250 
1.8636 
2 , il7Q6- 
5000 
0.^65217 
ff^f '^/^ p. 18958 

- 2.1933 



y ^0. 75000 
y = ^0.81818 , 

y = -i.im 

y -2.5000 . 

- y^f. -l'-2609 

' = -0.3397V 
= .0.30269 



J 



.Although this program will soive our problems, je have actually 
computed much more than neces^aiy. 5ie next section of J;he text 
shows how to solve equations ii^re efficiently. 



100 ' 



ERIC 
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Wl-h Simultaneous linear - equations : dev^mng a systemati c mejho^ ^ 

of solution 
Answers to Sjcercises F7-^ 



1. 



2. 



I READ 1/ AU,A12<. Bl, A21, A22,^B2 
' fRim 1. A1l/ai2, B1, A21, A22, B2 
FOR?;lATTip9-0^5FiO-0) ' 
AAli = An/All 
^AAl2.y^ A12/A11 . 
'SBl/^Bl/All 

fl = A21 - A21 ^ AAli 
/^2 .= A22 - AA12 • 

BS2 ^ B2 - A21 * BBl 
X2 = BB2/AA22 
Xi = BBl - AAX2 * X2 
PRINT 2, XI, X2 
2 - F0EMAT(THca3X7a=E20.8, 5Xl+HX2a= .E20.8) 

GO TO 3 



(a) -x = 1.6250 

(b) X*- 1.8636 

(c) ' X = 2.i;70>^ 
'(d) X =^".5000 

{ey<j'^ 0.^65217 
^th'" ,^.18958 

0\ - 2-1933 



y i"^0.75000 
y = ^0.81818 . 
y = -l.lkn 
y = -2.5000 . 
y_= -1:2609 
^ = -0.3397V 
= -0.30269 



J 



.Although this program vill solve our T)roblems, vj.e have actually 
computed much more than necessary. 5ie next section of J;he text 
shows how to solve equations mbre efficiently. ' , 



ICO ' 
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1F7-5- Simultaneous .linear equations : Gauss '- algorithm 
: ' Answers 'to Exercises Set A 

I 

1. 



Dd loe j'= 3/3 

: A(2,J) = A(2,JVa(2,2) 
100 - COmMJE 

B(2) = B(2)/A(2y2) 



DO 100 K-^ 1,2 
' . , . JX = K + 1 

- DO 200 J = J\,2 
A(K,J) = A(K,J)/a(K,K) 
200 *]COmNUE ■ * 
^ ^m) - B(K)/a(K,k; 
100 CONTINUE- 




' IF (K - 3) 5, '60, 60 • 
5. KPl = K + 1 ' 

DO 500 I = KPl, 3 • . ' 

• ; DO 600 J = KPl '3 

^A(I,J.) . A(I,J) . A(I,K)-* A{K,J) 

6eo coroiNUB . 

^'B(I)'= B(I) - A(I,K)'* B(K)' 
500 . CONTINUE • 



600 
^0 



IDO '100 K = 1,3 \ ' ' 
KPi = K + 1 ' . 

Da^2oo J = m,3. . . 7 

^ A(K,J) = A(K,j)/A(K,iC) 
- CONTINUE ' ' 
B(K) = B(K)/A(R,it) . ■ 
IF (KPl-liy?, ^0, ;60 V 
DO 500 I^= fcpi, 3 
DO 600 J = KPl, 3 

A(I,J) = A(I,J) - A(I,K) tUlK.JM^' 
^ CONTINUE / , - • * " ^ 

^I) = B(I) . A^I>4C) * b(k1 . . , , 
CONTINUE 



Answers to Exercises^F7^5 Set B 

1. 



i 

^2. 



31 
150 



200 



i 

6(k) 



1100 
900 



DIMENSION A(3,3), B(3)> X(3) 
DO 150'l'*^!l,3 

READ 3ir'\A(i;j)^ J = l,3)r .B(l) 
PRINT 31, '(A(I,J), J = 1;3)/ B(I^) 
FORMAT (iiFl'0;5) 
CONTINUE 
DO 100 K = 1X3 

KPl. = K -f 1" * * 

IF (K - 3) '^j-^, ^ 
DO 200 J = KL,^ 

A(K,j) = a(k:vJ)/a(k,k) 

CONTINUE \ 
B(K) = B(K)/A(K,K) 

IF (K - 3) 5/- 100, J.00 
DO 500" I = Ktl,3 • 
DO 600 J = KPr,3 \ 

A(l,h = A(I,J) V A(I,K) * A(K>J) 
W CONTINUE ' ' / * \ ' 

B(i/) =*B(i) - a(i\k)-* B(K) 

. CONTINUE • \ 

CONTINUE ,\ , 

DO 90b J = 1,3 ' \ ' ' 

. II = U -'l ; / \ ' . 
,X(I1) = B(I1) \ 

IF (I - 1) 11, "900, ii 

IlMl = 11-1 X 
DO 1100 J = 1, IlMl 'V * 
Jl = 1^ {I • • . f 

^•x(ii)r^ xiXD A(ii;4) * ^^"^^^^ 

CONTINUE ^ . •% . 

■ CONTINUE ' - '>'•'! ' '* ' ' 
PRINT 150, (X(l), I = i,3 \ . 

4 ^ 



STOP 
END 



-10 



.000,^ 



Xg = 1.882!+ 
= 15.'+71 



\ 

\c) = 1.2289 
\ x' = 0.20482 . 

i c. 

-Q. 83135 



2.6279 

-0'.2325'6 



x^ = -1.8372 



•(4) 



2.8154 

i'.TOTT' 
-0.15385 



-• (a), x^;= -3-0619 
> _ X2= 5.9268 
' \ = 0.13861 



>1) 



^1 



0.€63ll 

5. mi 

-1-5221 



Ansvers -to Exercise F7-$ Set C 



200 

4 

5 

6oo 
500 

100 



1100 
900 



SUBROUTINE GAUSS(N, A,B,xj 
DIMENSION A(15,15)., B(l5), X(l5) 
DO 100 K = 1,N I 
*KP1 .4- 1' 
IF (K . N) 2, h^k 
DO 200 J = Kl^N ^ 
' A(K,J) = A(K,J)/A(R,K) " ' 

CONTINUE 
B(K) = B(K)/A(K,K) ^ 

IF (K.. N). 5; 100, 100 

DO 500 I KPf;N- 

DO 600 J ^= Kri^N . 

A(I>J) =.A(r,j) - A(1,K) * A(K,J) 
CONTINUE 

•^9(1) - B(I} . A(I,K) ^ B(K) ' ' 
CONTinUE • • 

CONTINUE ' . 
DO 900^1 = 1,N - , ' • 
II = N -f 1 - I 

X(li) ='B(I1> . ' ^ 

IF (l-l) 11, 900; 11 

IlMl = Il-l 

DO 1100 J = 1, im 

Jil = N -f 1 - J 
X(ll) = X(I1) - A(I1,J1) ^ X(J1) 
CONTINUE . ' 

CONTINUE . ^ 

RETUiiN 



* I 

Here is a sample bailing ^rograiti. 



100 



101 
150 



DIMENSION A(15,15), B(l5)^ X(l5> 

RE^D 100, N t 

FORMAT (1$) . 

DO. 150 I = 1,N . ^ 

BEAD 101, (A(I,J), ^ = l.,N) 
PRINT 101, (ACl,J)^ J =:'1,1T) 
PORMAO? (5F10-5)' * • 

CONTINUE 

READ 101 (B(I), I = l^N) 
BRINT 101 (B(I), I,= 1,N). 
CALL GAUSS (N,A,B^) 
PRINT- 101, (X(I), I :=v^l,N> . ^ 
STOP . ' ' ^ > 
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Arisvers to gxercises F7-5 >^ Set B 
1. 



J 



(a) We add the index L as fifth argument f or indicating an alternate 
exit in' the ^ event the .'matrix A f is singul^V. i j^^^ 



(b)* 'Immediately after the statement ^ 

' ' % KPl + 1. 

we noV inseH the* following statements; 



li^OO 
22 

♦ • l8 



1900 



TOAX = ABSF(A(K,K)) 
M = K 

DO li+00 I = KP1,N ^ , * ' 
IF 'ABSF(A(I,K) - WX)ikOO,l6',l6 
MAX = ABSF(A(i;iC)) 
M = I 

CONTINUE • ' 

IF (MAX) 18,22,18 
I> = 1 * ' . 
RETUBN 

IF (M-K) 19>2A9» . .... 
DC 1^00 J W K,t • * • 
> COPY = A{K,J) . 

,A{K,J) =. A(M,J) • 
X(M,J) = COPY * . 

CONTINaE 
COPY = B(K) . 
^'BiKY = B(My 
b'(M) = COPY 



■2. Results vfith partial pixeting/ 



' (a) = 1.1805 

: X- .= 0.59398 

^ . 3 ^ 

Cbj/ = 10.550- 

■ - ' ■ X2 = 3.9000 

x^- = -0.60000 

• 3 ■ 



ERIC 



Without partial- pivpting 



(a-) = 



Xg = O.O92593 
1.7963 



Jb) x^ = 2.75' 
= 0 ■ . 
X3.=3.75 



Nate: Of course the ansyers without pivot are Wrong sincfe the machd^ne ^ 
continued tfte calcula^ after the attempted divisi6ii ^ zexo. Some^ 
compilers may perfojnn an error, stop if di>visioh by zero is encountered. 



4*' 



102 



